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FOREWORD 


Tnis  report  is  the  final  technical  report  on 
work  carried  out  under  Contract  DAHC04-70-C-0022 , 
which  was  administered  by  the  U.S.  Army  Research  Office, 
Durham.  The  findings  in  this  report  are  not  to  be 
construed  as  ar.  official  Department  of  the  Army  position 
unless  so  designated  by  other  authorized  documents. 


Abstract 


The  Carbon  Monoxide  laser  is  one  of  the  most  promising 
systems  currently  available  for  the  efficient  generation  of 
infrared  radiation  in  the  wavelength  region  from  5  to  8  urn. 

A  detailed  model  of  the  kinetic  processes  which  govern  the 
performance  of  this  laser  has  been  programmed  for  the  IBM 
360/67  computer.  The  temporal  development  of  vibrational 
populations  has  been  computed,  and  shown  to  converge  to  a  steady 
state  solution  generated  elsewhere  from  algebraic  computation. 

A  careful  investigation  of  scattering  data  shows  an  apparent 
discrepancy  in  the  potential  used  to  describe  the  interaction 
of  the  CO  molecule  with  Argon  atoms.  An  experimental  facility 
has  been  developed  for  observation  of  spontaneous  emission 
from  various  molecular  transitions  in  the  CO  laser. 
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I .  Introduction  and  Summary 

A  contract  was  awarded  to  Wayne  State  University  for  "A 
Study  of  Energy  Transfer  Processes  in  Molecular  Lasers",  with 
an  effective  starting  date  of  26  January,  1970.  The  effort 
under  this  contract  is  twofold.  A  computational  effort  is 
directed  at  modeling  energy  transfer  processes  important  to 
molecular  laser  action  using  detailed  computer  codes,  some  of 
which  were  already  in  operation  at  Wayne  State  University. 

The  experimental  part  of  the  effort,  which  is  relatively  small 
due  to  funding  limitations,  is  directed  at  experimental 
verification  of  the  computational  and  theoretical  assumptions 
implicit  in  the  calculational  model.  The  principal  investigator 
for  this  research  is  Dr.  Alexander  J.  Glass,  Professor  of 
Electrical  Engineering.  Other  faculty  associates  are  Drs. 
Richard  Marriott  and  Edward  Fisher,  Professors  of  Chemical 
Engineering.  All  three  are  members  of  the  Research  Institute 
for  Engineering  Sciences  at  Wayne  State  University. 

1 .  1  Description  of  the  Model 

Among  the  most  important  of  the  molecular  laser  systems 
currently  in  use  is  the  CO  laser.  This  laser  operates  between 
successive  vibrational  levels  of  the  CO  molecule,  which  is 
excited  in  an  electrical  discharge,  either  pulsed  or  CW,  in 
the  presence  of  He,  Nj ,  and  other  additive  gases,  most  usually 
02  and  Xe.  By  varying  the  conditions  of  excitation,  the  gas 
mixture,  and  the  dispersive  properties  of  the  laser  cavity, 
various  ( Av  =  1)  vibrational  transitions  can  be  stimulated, 
giving  the  CO  laser  a  wavelength  diversity  far  greater  than 
other  molecular  lasers.  As  examples  of  the  variation  in 
output  which  can  be  achieved  in  the  CO  laser,  the  range  of 
operating  wavelengths  which  have  been  observed  at  room 
temperature,  with  CW  oscillation,  by  Bhaumik  and  co-workers ^ ^ 
as  shown  in  Table  I  -  1.  Similarly,  in  Table  I  -  2,  we  show 
the  range  of  transitions  observed  with  pulsed  excitation. 
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at  77°K,  by  Graham  et  al.^)  __  it  is  the  main  purpose  of  the 
research  outlined  in  this  report,  to  develop  a  model  for  the 
CO  laser  system  which  is  of  sufficient  flexibility  to  provide 
a  sound  analytic  basis  for  determining  the  optimum  operating 
conditions  for  the  CO  laser,  for  either  pulsed  or  CW  operation, 
at  any  wavelength  in  the  range  of  operation. 

Stimulated  emission  in  any  laser  system  is  the  result 

( 3) 

of  the  creation  of  an  inverted  population  ,  a  condition  in 

which  the  population  of  a  given  excited  quantum  state  exceeds 

that  of  a  quantum  state  of  lower  energy  (barring  degeneracies) . 

In  order  to  achieve  inversion  in  the  CO  system,  several 

different  energy  transfer  processes  must  play  a  role.  To  begin 

with,  excitation  of  the  CO  molecules  by  impact  with  energetic 

electrons  in  the  discharge  must  take  place.  As  is  shown  in 

Chapter  III  of  this  report,  there  are  no  selective  excitation 

processes  taking  place  in  the  CO  laser,  at  least  not  to  any 

great  degree,  so  that  from  excitation  rates  only,  we  would 

not  expect  to  find  an  inversion  of  population.  It  is  not  until 

we  consider  the  effect  of  vibrational  energy  transfer  (V  -  V) 

that  the  inversion  mechanism  becomes  clear.  In  an  an  harmonic 

system  like  CO,  providing  that  the  average  kinetic  energy  of 

the  molecules  is  kept  well  below  the  vibrational  level  separation, 

a  significant  deviation  from  a  Boltzmann  distribution  of 

vibrational  populations  can  be  created,  with  an  accumulation  of 

(4) 

population  in  the  higher  lying  states  .  Detailed  calculations' 
of  this  effect  are  exhibited  in  Chapter  III  of  this  report. 

j 

In  order  to  assure  that  the  gas  temperature  remains  low,  helium 

J  (  I 

is  added  to  the  discharge,  to  provide  good  heat  conduction  tp 
the  walls.  Thus  the  model  must  include  electron  excitation  of 
vibrational  states,  V  -  V  and  V  -  T  processes  among  the  molecular 
species  present  in  the  laser,  and  the  effect  of  neutral  atom 
collisions  with  the  molecules.  Each  of  these  processes  is 
characterized  by  a  cross  section,  which  depends  on  the  relative. 

,  i 
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velocity  of  the  collision  partners.  These  cross  sections 

must  be  integrated  over  the  appropriate  velocity  distributions 

to  obtain  the  rate  coefficients  for  the  respective  processes. 

It  is  assumed  in  this  research  that  the  kinetic  energy 

distributions  of  the  atomic  and  molecular  constituents  of 

the  gas  are  characterized  by  a  single  gas  temperature,  and 

are  described  by  a  Maxwell-Boltzmann  distribution  of  velocities. 

The  vibrational  levels  are  treated  explicitly,  by  a  set  of 

coupled  rate  equations  for  the  vibrational  populations,  the 

master  equation.  The  electronic  velocity  distributions  are 

assumed  to  be  Maxwellian,  at  a  temperature  different  from 

( 5 ) 

the  gas  temperature.  It  has  been  shown  that  in  a  weakly 
ionized,  molecular  plasma,  the  electron  velocity  distribution 
deviates  from  Maxwellian  form,  but,  as  is  discussed  in 
Chapter  III,  this  deviation  has  little  effect  on  the  vibrational 
excitation  rate  coefficients. 

Radiation  processes  must  also  be  considered  in  the 
model.  Only  spontaneous  emission  has  been  treated  to  date, 
so  the  model  only  applies  to  laser  performance  below  and  up  to 
the  threshold  of  oscillation.  This  is  sufficient  for  comparison 
with  observation,  however,  since  the  principal  diagnostic 
measurement  is  spontaneous  emission  from  the  various  levels 
in  the  discharge  (sidelight  emission) . 

I .  2  Summary  of  the  Report  and  Results 

Chapter  II  of  the  report  considers  the  vibrational 
excitation  of  carbon  monoxide  by  collisions  with  inert  gas 
atoms,  in  particular,  argon  and  helium.  In  both  cases, 
experimental  data  on  vibrational  relaxation  of  CO  in  the 
presence  of  the  atomic  gas  are  available,  for  comparison  to 
the  calculated  values.  A  detailed  description  is  given  of 
the  close  coupling  calculation.  It  is  necessary  in  this 
calculation  to  represent  the  atom-molecule  interaction  by 
sore  potential  function  of  the  separation  coordinate.  At 

/  g\ 

first,  the  Lennard-Jones  (6  -  12)  potential'  was  employed. 
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The  computed  relaxation  times  were  smaller,  by  significant 

factors,  than  experimental  data,  both  for  He-CO,  Ne-CO  and 

Ar-CO.  In  order  to  determine  the  cause  of  this  discrepancy, 

the  calculations  were  repeated  using  an  empirical,  "soft-core" 

( 7) 

potential  ,  obtained  from  scattering  data.  Using  the 
modified  potentials,  satisfactory  agreement  with  experiment 
was  obtained  for  CO-He  mixtures,  but  for  CO-Ar  mixtures,  the 
calculated  relaxation  times  exceeded  the  experimental  values 
by  a  factor  of  50.  This  unexpected  disagreement  may  be  due 
to  the  presence  of  quenching  impurities,  such  as  1^0,  or  H2> 
in  the  experiment,  but  it  is  felt  that  this  is  unlikely.  The 
exact  cause  of  the  discrepancy  remains  unresolved  at  present. 

In  Chapter  III,  a  detailed  summary  is  given  of  the 

computational  approach  employed  in  the  solution  of  the  master 

equation.  Up  to  30  vibrational  levels  of  both  N2  and  CO,  or 

indeed  any  pair  of  diatomic  gases,  can  be  explicitly  included, 

along  with  electrons  and  inert  atoms.  A  preliminary  test  of 

the  code  was  performed  by  considering  the  reaction  of  vibrationally 

excited  H2  and  D2  to  form  HD  in  the  absence  of  electrons. 

Detailed  comparisons  are  given  for  the  Raman  excitation  of  the 

vibrational  levels  of  H0  in  the  presence  of  D0,  as  carried  out 

( 8 )  ^  ^ 

by  Bauer  et  al.  .  Good  agreement  with  experiment  is  obtained. 

It  is  also  pointed  out  that  a  more  definitive  reaction  in  order 
to  test  the  proposed  HD  mechanism  would  proceed  from  the 
excitation  of  D2  in  the  presence  of  H2.  To  date,  the  latter 
experiment  has  not  been  carried  out. 

The  dependence  of  the  vibrational  excitation  rates  on 

electron  temperature  is  then  discussed.  It  is  found  that  the 

excitation  rates  are  insensitive  to  the  shape  of  the  high  energy 

tail  of  the  electron  distribution  T  5s  10000°K  for  vibrational 

e 

excitation.  Accordingly,  a  Maxwellian  distribution  can  be  used 

to  describe  the  electron  distribution  in  the  model.  Finally, 

the  transient  response  of  the  CO-N2  laser  mixture  to  an  initially 

established  electron  temperature  is  computed,  and  is  shown  to 

( 9) 

converge  in  time  to  the  steady- state  results  of  Rich  . 
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In  Chapter  IV,  a  brief  di scription  is  given  of  the 
experimental  facility.  There  are  three  appendices  to  the 
report,  a  discussion  of  the  calculation  of  V  -  V  rate 
coefficients,  a  summary  of  the  modeling  calculations,  which 
was  presented  at  the  Conference  on  Laser  Physics  at  Sun  Valley, 
Idaho,  on  March  1-3,  1971,  by  Drs.  Fisher  and  Kummler,  and 
a  discussion  of  the  vibrational  excitation  of  CO  by  collision 
with  He,  a  paper  to  be  given  at  VII  International  Conference 
on  the  Physics  of  Electronic  and  Atomic  Collisions,  by 
Dr.  Marriott. 
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TABLE  I 

-  1 

Spectral 

Output  of  CO 

Laser  at  20°C(^m) 

X  Observed 

X  Theoretical  .. 

Vibrational 

Relative 

(Air) 

(Air) 

v  cm- 

Band 

Transition 

Intensity 

5. 3155 

5.3152 

1880.897 

8-7 

P  ( 20) 

80 

5.3275 

5.3273 

1876.629 

P  (21) 

9 

5.3397 

5.3395 

1872.329 

P  (22) 

10 

5.3871 

5.3876 

1855.615 

9-8 

P  (20) 

100 

*  5.3993 

5.3999 

1851.382 

P  ( 21) 

16 

l 

1842.821 

10-9 

P(17) 

5.4487 

5.4494 

1834.577 

P  (19) 

52 

5.5131 

5.5127 

1813.514 

11  -  10 

P  ( 18) 

17 

5.5256 

5.5252 

1089.416 

P  ( 19) 

1.8 

5.5901 

5.5901 

1788.398 

12  -  11 

P  (18) 

2.2 

5.6027 

5.6028 

1784.334 

P  (19) 

0.9 

5.6559 

5.6567 

1767.359 

13  -  12 

P  (17) 

0.3 

5.6688 

5.6695 

1763.363 

P(18) 

0.5  1 

5.6819 

5.6825 

1759.334 

P  (19) 

— 

5.6991 

5.6996 

1754.060 

14  -  13 

P  (14) 

0.3 

1 
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TABLE  1-2.  Partial  Summary  of  Observations 
at  77 °K  for  0.3  Torr  CO  Plus  5  Torr  Helium. 


Observed 

Transitions 

Strongest  Line  Data 

V  -  V' 

P(J) 

(a) 

B 

\  vac 
(microns) 

Current 

(Amp) 

Laser 

Output 

(watts) 

Laser 

Pulse 

(msec) 

4-3 

5  - 

8 

■  I 

4.90882 

0.198 

1.2 

0.2 

5-4 

5  - 

8 

7 

4.97283 

0.193 

0.68 

0.3 

6-5 

,  5  - 

8 

7 

5.03735 

0.193 

0.40 

0.3 

7-6 

i  5  - 

8 

8 

5.11437 

0.335 

0.47 

0.4 

00 

!  5  - 

8 

7 

5.17189 

0.285 

0.12 

0.1 

9-8 

6 

6 

5.23190 

0.213 

0.005 

0.3 

10-9 

i  7 

7 

5.31342 

0.713 

»  0.007 

5.0 

(a) 


Includes 


only  those  transitions  for  which  laser  power 


exceeded  1  mw. 
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II.  A  THEORETICAL  STUDY  OF  THE  VIBRATIONAL  EXCITATION 
OF  CARBON  MONOXIDE  BY  INERT  GAS  COLLISIONS 

1.  introduction 

Inelastic  processes  involving  the  vibrational  excitation 
of  molecules  by  collision  with  inert  gases  are  of  particular 
relevance  to  the  study  of  gas  laser  phenomena.  Much  of  the 
necessary  data  is  not  available.  Direct  experimental  measure¬ 
ments  of  the  corresponding  low  energy  collision  cross  sections 
have  not  been  made,  and  the  non  Boltzmann  nature  of  the  system 
may  well  render  invalid  any  extrapolation  based  on  laboratory 
observation  of  vibrational  relaxation  processes.  It  is  of 
importance,  therefore,  to  establish  the  quantitative  reliability 
of  the  theoretical  procedures  for  calculating  the  rate  processes 
required. 

In  previous  publications  (Marriott  ^  Gianturco  and 

/  O  \ 

Marriott  v  )  a  method  for  the  close  coupling  calculation  of 
cross  sections  for  the  vibrational  excitation  of  molecules  by 
collision  with  atomic  particles  has  been  developed  and  applied 
to  a  variety  of  collision  processes. 

In  this  report  the  method  is  applied  to  the  vibrational 

excitation  of  CO  by  collision  with  Ar  and  He.  There  are  no 

anomalous  features  in  the  interactions  between  the  particles 

in  either  case,  consequently  the  results  obtained  would  be 

expected  to  be  as  reliable  as  those  obtained  in  the  original 

(3) 

study  of  CO  -  CO  collisions  (Gianturco  and  Marriott 
Experimental  data  is  available  on  the  vibrational  relaxation 

(A)  (5) 

of  the  mixtures  (Hooker v  , Millikan  ')  for  comparison  with 
the  calculations. 
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As  summarized  in  section  2,  the  close  coupling  calculation 

is  based  upon  the  assumption  of  a  spherically  symmetric  semi- 

empirical  scattering  potential.  However,  it  is  found  that  use 

of  the  appropriate  Lennard-Jones  potential  as  was  done  in  the 

original  calculations,  leads  to  vibrational  relaxation  times 

that  are  too  short,  hence  by  inference  to  cross  sections  for 

the  excitation  of  the  first  vibrational  level  of  CO  that  are 

too  large.  The  details  of  these  calculations,  including  the 

result  for  a  similar  CO-Ne  calculation,  are  presented  in  sections 

3,  4  and  5.  Sections  6  and  7  present  a  recalculation  of  the 

cross  sections  for  excitation  by  Ar  and  He  based  upon  the  more 

(6) 

realistic  short  range  potentials  measured  by  Amdur  et  al . 
and  Jordan  et  al.  ^  for  collision  energies  up  to  2.5  eV. 

This  leads  to  satisfactory  agreement  with  experimental  relax¬ 
ation  times  for  the  case  of  CO-He  mixtures.  However,  for  the 
CO-Ar  system  it  is  found  that  the  calculated  relaxation  times 
are  now  greater  than  experiment  by  about  a  factor  of  50. 

This  surprising  disagreement,  occuring  in  what  would  have 
been  expected  to  be  a  routine  calculation,  is  investigated  in 
detail  in  section  8.  The  theoretical  results  are  confirmed. 

At  the  same  time  it  is  estimated  that  the  degree  of  contamination 
of  Ar  by  ^  or  ^0  which  could  account  for  an  experimental 
error  of  this  size,  though  small  would  still  be  two  orders  of 
magnitude  larger  than  is  normal.  The  disagreement  remains 
unresolved . 
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2.  Collision  theory 


In  summary,  the  analysis  of  the  collision  process  neglects 
the  effects  of  the  rotational  states  of  the  target  and  the 
internal  structure  of  the  impacting  particle.  It  assumes  the 
molecular  wave  function  to  be  separable  into  electronic, 
rotational  and  vibrational  components  and  the  scattering  potential 
to  be  similarly  separable  in  terms  of  the  molecular  normal 
coordinates  and  to  be  spherically  symmetric. 

To  this  approximation  the  coupled  radial  scattering 
equations  take  the  form  using  atomic  units  throughout. 


+  k 


2  _  iU+1) 
n  2 


r 


(1> 


where  the  wave  number  of  the  scattered  particle  of  reduced 
mass  M  is  given  by 

k2  =  2M(E  -  E  +  k2/2M)  (2) 

n  p  n  p 


being  the  energy  of  the  nth  molecular  vibrational  state  and 
the  initial  state  being  arbitrarily  assigned  the  subscript  p. 


The  scattering  equations  (1)  are  solved  numerically  with 
the  retention  of  all  coupling  between  any  specified  set  of 
molecular  states  subject  to  the  boundary  conditions 

F*(0)  =  0 
n 

F*  (r)~iA  (2  4-i  1 )  ["  6  k_1sin(k  r-t.Tr/2)  +  exp(ik  r-ilv/2)ft  1 
n  L  pn  p  p  n  pnJ 


(3) 
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The  corresponding  inelastic  partial  cross  sections  for 
the  collision  induced  transitions  are  then  given  by 


k 

4tt  ,-n.  (2R1) 

k 

P 


(4) 


and  the  total  cross  sections  obtained  by  graphical  integration 
over  all  significant  values  of  l. 


The  transition  matrix  elements  for  a  diatomic  molecule 
are  shown  to  be  approximated  by 


Ml  (r)  =  v  V(r )Vl 
nm  o  nm 


(5) 


where  numerical  values  for  the  factors  are  estimated  in  the 

nm 

usual  way  from  the  classical  distance  of  closest  approach  of 

( 2) 

the  colliding  particles  (Marriott  )  . 

Empirical  functions  are  employed  to  represent  the  scattering 

potential  component,  V(r),  and  the  normalising  constant  Vn  is 

o-l  u 

taken  to  be  (vqq)  .  This  ensures  that  the  scattering  potential 

for  the  ground  state  molecule  is  given  by  the  observed  empirical 

function . 

This  close  coupling  analysis  of  the  collision  process  has 
been  programmed  in  Fortran  for  the  IBM  360  series  computer. 


II 


4 


3.  The  Lennard-Jones  intermolecular  potential 


Previous  calculations  for  collision  systems  involving  CO 
have  successfully  employed  the  semi-empirical  Lennard-Jones 
functions  to  represent  the  scattering  potential  between  the 

(3) 

CO  and  the  impacting  particle  (Gianturco  and  Marriott  '  )  • 

To  this  approximation  the  potential  is  given  by  a  function 
of  the  form 


V(r)  --=  (f)6] 


where  the  necessary  parameters  can  be  estimated  using  the 
combin 
et  al. 


combination  rules  for  non  polar  mixtures  discussed  by  Hirschfelder 
(8) 


i. 

-i2 


c"«  =  (€coex  ’  '  ”  %te,co+  V1 

Hirschfelder  et  al.  give  for  the  component  parameters: 


Eco/k  = 

0 

CO 

CO 

°co  “  3-706 

X 

,  -8 

10  cm 

eArA  “ 

116°K 

°Ar  ^  3 • 465 

X 

,  -8 

10  cm 

eM  /k  = 
Nc' 

27 . 5°K 

°Ne  -•  2.858 

X 

10  ®cm 

l,  /k  - 
Ho 

10 . 2°K 

°fic  “  J-55f> 

X 

10  ^cin 

k  being  the  Boltzmann  constant,  and  these  lead  in  turn  to  the 
values  for  the  mixtures; 
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I 


CO 

H  Ar  : 

€  ~ 

.320  x 

] U~3a.u. , 

0  ~ 

CO 

H  Ne  : 

e  ~ 

.156  x 

10  3a . u.  , 

o  - 

CO 

•h  Ho  : 

L  ~ 

.950  x 

.  A 

10  'a  .  u .  , 

a  - 

0' 

V 

'0] 


A  typical  matrix  of  transition  elements  obtained  for  the 
CO  +  Ar  system  at  several  collision  energies  and  an  impact 

i  I 

parameter  of  zero  are  shown  in  table  1.  , 


i  i 


■i '  i 
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i 
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Table  1.  Matrix  elements  for  head-on  collision  of  Ar  with  CO 
assuming  a  Lennard-Jones  interaction. 


Collision 


Transition 

elements* 


!eV)  + 

0.6 

1.0 

1.5 

2.0 

2.5 

o 

o 

> 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

V10 

0.0777 

0.0801 

0.0820 

0.0834 

0.0846 

V11 

1.0060 

1.0064 

1.0067 

1.0070 

1.0072 

V20 

0.0043 

0.0045 

0.0048 

0.0049 

0.0051 

V21 

0.1102 

0.1137 

0.1164 

0.1184 

0.1200 

V22 

1.0121 

1.0129 

1.0135 

1.0140 

1.0143 

V30 

0.0002 

0.0002 

0.0002 

0.0002 

0.0002 

V31 

0.0074 

0.0079 

0.0083 

0.0085 

0.0088 

V32 

0.1345 

0.1396 

0.1430 

0.1455 

0.1475 

V33 

1.0182 

1.0193 

1.0203 

1.0210 

1.0215 

+Relative  to  ground  state 


*VQ0  nomalized  to  1.0 
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4.  Vibrational  excitation  cross  sections  (LJ) 


Using  the  Lennard- Jones  potential  described  in  the  previous 
section  the  remaining  parameters  required  as  data  for  the 

5 

computer  code  are  the  reduced  oscillator  mass  (0.126  x  10  atomic 
units),  the  vibrator  energy  level  separation  for  CO  (0.266  eV) 
and  the  collision  reduced  masses: 

5  5 

CO  +  Ar;  .3024  x  10  a.u.,  CO  +  Ne;  .2149  x  10  a.u., 

4 

CO  +  He;  .6433  x  10  a.u. 

The  calculated  partial  sections  for  the  excitation  of  the 
first  vibrational  state  by  Ar  are  shown  in  figure  1  where 
coupling  between  the  first  four  states  has  been  retained  and 
in  figure  2  where  coupling  has  been  restricted  to  the  first 
two  states.  The  total  cross  sections,  obtained  by  graphical 
integration  of  these  figures  are  presented  as  functions  of 
collision  energy  in  table  2. 

An  indication  of  the  numerical  accuracy  attained  was  given 
by  the  fact  that  detailed  balance  held  in  all  cases  to  at  least 
3  significant  figures.  Inclusion  of  more  than  four  vibrational 
states  was  not  found  to  have  any  further  significant  effect  on 
the  Qq1  cross  section. 

The  results  for  the  Ne  and  He  calculations  were  quite 
similar  in  form  to  those  for  Ar  and  the  details  are  not 
included  here. 

Cross  sections  for  the  excitation  of  CO  from  the  ground  to 
the  first  vibrational  level  in  the  three  collision  systems  are 
shown  in  table  3.  In  these  calculations  coupling  has  been 
retained  between  the  first  3  or  4  vibrational  state  of  CO. 
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PARTIAL  WAVE ,  £ 

Figure  1.  Calculated  partial  cross  sections  for  the 
excitation  of  the  first  vibrational  level  of  CO  by 
collision  with  Ar  for  a  Lennard— Jones  interaction  at 
several  collision  energies  and  including  coupling 
between  the  first  four  vibrational  states. 

A,  collision  energy  0.8  eV;  E,  collision  energy  1.0  eV; 
C,  collision  energy  1.5  eV.. 
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(ndy 


0  40  60  120  160  200  240 

PARTIAL  WAVE,  t 


Figure  2.  Calculated  partial  cross  sections  for  the 
excitation  of  the  first  vibrational  level  of  CO  by 
collision  with  Ar  for  a  Lennard- Jones  interaction  at 
several  collision  energies  and  including  coupling 
between  the  first  two  molecular  vibrational  states 
A,  collision  energy  0.6  eV;  B,  collision  energy  1.0  eV; 
C,  collision  energy  1.5  eV?  D,  collision  energy  2.0  eV; 
E,  collision  energy  2.5  eV. 
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Table  2.  Total  cross  sections  for  vibrational  excitation  of 

CO  by  Ar  impact  assuming  a  Lennard-Jones  interaction 
potential. 


Collision 

Two  coupled  states 

Four  coupled 

energy  (eV) 

°01(’,“o) 

Q01 ( nao 

0.266 

0.00 

„  -7* 

0.00 

0.6 

0.20 

-5 

0.8 

— 

0.98 

1.0 

0.35"4 

0 . 16~3 

1.5 

0 . 20-2 
_2 

0.57“2 

2.0 

0.14 

- - 

2.5 

0.42-1 

_  _ _ 

Table  3.  Total  cross  sections  (LJ)  for  the  vibrational  excitation 
of  CO  from  the  ground  to  the  first  excited  level  by 
collision  with  Ax",  Ne  and  He. 


Collision  energy  (eV) 

0.6 

0.8 

1.0 

1.5 

to 

• 

o 

Helium 

Q01  (nao) 

. 16-1* 

• 

00 

•p* 

.24 

.97 

.19+1 

Neon 

Q0l<Tra2) 

.24"4 

.  55-3 

.39-2 

.  66-1 

.30 

Argon 

- 

-5 
.98  3 

.  16”3 

.  57-2 

- 

*  superscript  denotes  the  power  of  10  by  which  the  number  must 
be  multiplied. 
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5.  Vibrational  relaxation  of  CO  in  the  inert  gases 


No  direct  measurements  of  vibrational  excitation  cross 
sections  for  heavy  particle  collisions  have  yet  been  made. 
Comparison  of  theory  with  experiment  can  only  be  obtained  in 
an  indirect  way  by  using  the  cross  sections  to  determine  the 
relaxation  time  for  CO  in  the  corresponding  gas.  These  relax¬ 
ation  times  may  then  be  compared  with  the  experimental  results 
of  Hooker  and  Millikan  ^  ^  and  Millikan ^  . 

The  connection  between  the  relaxation  time,  t,  for  a 

simple  harmonic  vibrator  and  the  corresponding  process  reaction 

rate  coefficient,  v  »  has  been  discussed  by  Herzfeld  and 
(9)  nm 

Litovitz 

If  it  is  assumed  that  vibrational  transitions  can  occur 
only  to  neighboring  states,  then  it  can  be  shown  that,  for  a 
diatomic  gas  in  translational  equilibrium 


1/t  =  v10  "  V01 


(6) 


where  in  general  for  a  gas  mixture  of  C  components  of  fractional 
molecular  concentrations 


(7) 


where 


4ttN 


(v)cxp(-  M^vs/2kT)  v3dv 


N  being  the  number  density  of  the  gas,  g^  the  statistical 

weight  of  the  m  vibrational  state,  M^  the  reduced  mass  for  the 

molecular  vibrator  and  the  i ' th  mixture  component  and  the 

Q1  (v)  the  cross  section  for  the  n-»m  transition  induced  by 
nm 

collision  at  velocity  v  with  the  i ' th  mixture  component.  For 


II  -  12 


the  reverse  reaction,  detailed  balancing  gives  immediately 


v  g 

mn  m 


exp  ;-f:  At)  =  v  q 
m  nm  n 


exo  (-E  At)  . 
n 


The  relative  magnitudes  of  the  transition  elements  in 
table  1  indicate  that,  as  was  the  case  for  pure  CO,  the  assump¬ 
tion  that  only  single  quantum  jump  transitions  can  occur 
should  be  valid  to  a  good  approximation  for  CO  in  Ar.  A 
similar  conclusion  is  reached  for  the  other  inert  gases. 

The  relaxation  time  for  CO  in  Ar  obtained  on  this  basis 
by  the  numerical  integration  of  equation  (7)  using  the  excit¬ 
ation  cross  sections  in  table  3  and  taking  N  to  be  the  particle 
density  for  1  atmosphere  pressure  at  T°K,  is  shown  as  a 
function  of  temperature  in  figure  3.  Substitution  of  the 
parametric  values  for  the  CO-Ar  system  reduces  equation  (7)  to 

/Qnm  (E)  exP  (”eAt)  (EAT)d(EAT)  sec  1  atn'“ 

(8) 

.  2 
where  Q  is  in  units  of  ira  . 
nm  o 

The  experimental  results  of  Hooker  and  Millikan  are  also 
shown  in  figure  3.  It  is  immediately  apparent  that  the  theor¬ 
etical  relaxation  time  is  consistently  smaller  than  the  observed 
value  by  a  factor  of  four.  This  in  turn  implies  that  the  calc¬ 
ulated  cross  sections  for  the  excitation  of  the  first  vibrational 
level  of  CO  by  collision  with  Ar  are  too  large  by  the  same 
factor. 

For  the  relaxation  of  CO  in  He  the  required  rate  coefficient, 
in  the  same  notation,  is  given  by 

0  GU 

Y0r=  4(10)n(MT)"Mo03  (E)  (E/kT)e“E/kTd(K/kT)  sec"1  atm. 

(9) 


the  convenient  form 

2.303  x  109 

Ynm  & 

T* 
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Figure  3.  Vibrational  relaxation  time  for  CO  in  Ar 
"at  one  atmosphere  pressure. 

(4) 

A.  Experimental  results  of  Hooker  and  Millikan 

B.  Calculated  in  this  paper  using  a  Lennard-Jones 

scattering  potential  with  coupling  between  the  first 

% 

four  vibrational  states. 


The  theoretical  relaxation  times  for  CO  in  He  obtained 
in  this  way  are  shown  in  figure  4  in  comparison  with  the 
experimental  values.  It  is  again  apparent  that  theory  is 
consistently  less  than  experiment,  now  by  a  factor  of  fifty. 

A  similar  result  is  obtained  for  CO  in  Ne  which  is 
intermediate  between  those  for  Ar  and  He.  The  theoretical 
relaxation  time  for  the  Ne  mixture  is  found  to  be  low  by  a 
factor  of  seven. 

In  view  of  the  excellent  agreement  between  theory  and 
experiment  for  the  case  of  pure  CO  (Gianturco  and  Marriott 
these  discrepancies  appear  to  be  significant. 

It  is  of  importance,  therefore,  to  establish  where  the 
approximations  that  proved  satisfactory  in  the  treatment  of 
CO-CO  collisions  fail  for  the  similarly  weakly  interacting 
inert  gas  systems.  Evidently  the  part  of  the  present 
calculation  that  involves  greater  approximation  than  was  the 
case  for  pure  CO  is  the  specification  of  the  scattering 
potential  itself. 

Taken  together  the  results  indicate  that  the  repulsive 
core  of  the  Lennard-Jones  potential,  which  dominates  the 
vibrational  excitation  process,  seriously  overestimates 
the  strength  of  the  short  range  interaction  between  CO  and 
the  inert  gases. 
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Figure  4.  Vibrational  relaxation  time  for  CO  in  He  at 
one  atmosphere  pressure.  Experimental  results  of 
Millikan  #  Theory  (LJ) ;  calculated  in  this  paper 

using  a  Lennard-Jones  potential  with  coupling  between 
up  to  4  vibrational  states. 
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6.  The  soft  core  intermolecular  potential 

/  C\ 

The  scattering  studies  by  Amdur  et  al.  v  and 
( 7) 

Jordan  et  al .  of  a  number  of  gas  molecule  potentials 

indicate  that  the  repulsive  part  of  the  interaction  between 

-7  -12 

CO  and  the  inert  gases  behaves  as  r  rather  than  r 

This  represents  a  much  softer  interaction  than  the  Lennard- 

_  12 

Jones  core,  which  behaves  as  r  ,  and  thus  it  follows  that 
use  of  the  Lennard-Jones  potential  will  yield  cross  sections 
that  are  too  large. 

For  the  CO-Ar  interaction  the  soft  core  potential  is 
given  in  the  form 

V(r)  =  47529/r6* 99  (10) 

for  particle  separation  in  the  range 
3.8  *  r(aQ)  £  5. 

Since  no  experimental  information  is  available  for  , 
separations  greater  than  this  it  is  assumed  in  the  following 
calculations  that  the  Lennard-Jones  potential  discussed  ip 
Section  3  can  be  used  to  represent  both  the  long  range  inter¬ 
action  outside  the  potential  minimum  and  also  the  location 
of  the  potential  zero. 

Subject  to  these  restrictions  a  polynomial  fit  was 

•  }  ' 

constructed  which  smoothly  joined  the  short  range  soft  core 
to  the  long  range  Lennard-Jones  function.  •  i 

The  specific  analytic  expression  used  was  . 

i 

V(r)  =  . 0422 ( a  -  r)  +  .0684(a  -  r) 2  +  .0422(a  -  r) 3 

+  .00692  {n  -  r)4  (eV)  ,  (  ('ll) 

where  i 

5  *  r(aQ)  s  7.5 

and  n  is  the  appropriate  Lennard-Jones  parameter. 

I 

1 
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For  separations  r(a^)  >  7>  5  the  Lennard-Jones  function 
described  in  Section  2  was  (used. 

'  I  ,  ?  } 

Figure  5  shows  schematically  the  complete  function. 

Clearly  this  potential  cannot  be  validly  employed  for  collision 
energies  greater  than  4.2  eV.  However,  this  is  considerably 
higher  than  the  energies  of  interest  here. 

1  ,  i 

1  i  ■ 

The  matrix  of  transition  elements  obtained  on  the  basis 
of  this  potential  at  several  Collision  energies  and  zero  impact 

I  I  ■  ! 

parameter  are  shown  in  table  4 .  ! 

Comparison  of  table  4  with  table  1  shows  immediately  that 
for  the  soft  core  potential,/  as  for  the  LSnnard- Jones ,  transitions 
from  the  ground  to  the  first  excited  vibrational  level  would  be 
expected  to  dominate  the  vibrational  relaxation  procesls.  In 
addition  corresponding  matrix  elements  are  now  smaller  by 
a  factor  of  about  2.  ( 

I  1  !  ! 

For  the  CO-He  system  the  interaction  has  not  been  measured 
directly.  However,  data  is  given  on  the  N„-  He  system  ,  and  it 

,  '  !  .  ,  !  .  I 

is  shown  that  the  CO  and  N~  interactions  are  very  similar. 

'  ,  Z  •  ■  I  !  '  I 

On  this  basis  we  estimate  the  repulsive  interaction  for 

i  n 

CO-fJe  to  be  232/r  (a.u.)  for  separations  3  <  r(aQ)  <4.3. 

/For  tlr^e  long  range  part  the  Lennard-Jones  potential  is  assumed 
still  valid  and  the  two  functions  are  smoothly  joined  over  a 
small  region  by  the  commone  tangent:  1  , 

/  1  ,  i  i 

V(r)  ,=  .04739  -  .  009158r  for  4.53  <  r(a  )  <  4.84. 

,  I  1  C  :  1  I 

Use  of  this  soft  core  (SC)  potential  again  reduces  the 

,  ! 

strength  of  the  coupling  between  vibrational  states  by  about 

i  1  '  . 

a  factor  of  2  as  was  the  case  for  the  CO-Ar  soft  core  potential. 

•  *  1  1 

This  reduction  in  magnitude  of  the  cdupling  terms  would  be 
expected  to  lead  in  turn  to  a  reduction  in  rhagnitude  of  the 
cross  sections.  1  ’  1  , 

i  ,  i 


/  / 
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V(r)  (eV) 


Figure  5.  Empirical  soft  core  CO-Ar  interaction  potential 


Region 

Region 

Region 

Region 


t i)  Potential  undefined 

i 

(ii)  Empirical  potential  after  Jordan  et  al. 
fiii)  Polynomial  fit 

(iv)  Lennard- Jones  after  Hirschfelder  et  al. 


Table  4  .  Transition  elements  for  head  on  collision  of 
Ar  with  CO  and  the  soft  core  interaction. 


Collis ion 

Energy  (eV) 

0.6 

1.0 

1.5 

2.0 

2.5 

o 

o 

> 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

V10 

.0421 

.0460 

.0493 

.0517 

.0535 

V11 

1.0018 

1.0021 

1.0024 

1.0027 

1.0029 

Transition  V2o 

.0013 

.0015 

.0017 

.0019 

.0020 

elements* 

V21 

.0596 

.0651 

.0698 

.0732 

.0758 

V22 

1.0035 

1.0042 

1.0049 

1.0057 

1.0057 

V30 

.0000 

.0000 

.0000 

.0001 

.0001 

V31 

.0022 

.0026 

.0030 

.0033 

.0035 

V32 

.0730 

.0798 

.0856 

.0897 

.0930 

V33  1.0053  1.0064 

+  Relative  to  ground  state 
*  VQ0  normalized  to  1.0 

1.0073 

1.0080 

1.0086 
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7.  Excitation  cross  sections  (SC)  for  CO  collision  with  Ar  and  He 

The  partial  cross  sections  for  CO-Ar  collisions  calculated 
on  the  basis  of  the  corresponding  soft  core  potential  described 
in  the  previous  section  and  the  remaining  collision  parameters 
listed  in  section  3  are  shown  in  figures  6  and  7.  The  total 
cross  sections  obtained  by  graphical  integration  of  the  curves 
for  the  case  of  two  and  three  coupled  vibrational  states  are 
shown  in  table  5  together  with  the  corresponding  Lennard-Jones 
cross  sections  for  excitation  from  the  ground  to  the  first 
vibrational  level  for  comparison. 

As  expected  the  (SC)  cross  sections  are  smaller  than  their 
corresponding  counterparts  by  about  a  factor  of  2  at  low  energies. 
However,  this  difference  increases  with  energy,  being  about  an 
order  of  magnitude  at  1  eV  and  rising  to  a  factor  of  50  at  2.5  eV. 

The  divergence  between  the  (SC)  and  (LJ)  cross  sections 
increases  with  collision  energy  as  a  consequence  of  the  rapid 
increase  of  the  difference  between  the  potential  functions  as 
the  classical  distance  of  closest  approach  decreases. 

The  results  of  the  CO-He  calculation  were  quite  similar 
in  form  to  those  for  Ar  although  larger  by  orders  of  magnitude, 
and  are  not  presented  in  detail  here.  The  total  cross  sections 
obtained  for  the  CO-He  system  are  shown  in  table  6,  once  again 
with  the  corresponding  Lennard-Jones  results  for  comparison. 

The  numerical  convergence  of  these  He  cross  sections  was  confirmed 
by  retaining  coupling  between  up  to  six  vibrational  states. 

Overall  numerical  accuracy  was  indicated  as  previously  by  the 
fact  that  detailed  balance  was  satisfied  throughout  to  3 
significant  figures. 
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PARTIAL  WAVE , f 

Figure  6.  Partial  cross  sections  for  the  excitation  of 
CO  from  the  ground  to  the  first  vibrational  level  by 
Ar  impact,  calculated  at  several  collision  energies  •' 
assuming  a  soft  core  scattering  potential  and  with 
coupling  retained  between  the  first  two  vibrational 
states.  A,  0.6  eV;  B,  1.0  eV;  C,  1.5  eV;  D,  210  eV; 

E,  2.5  eV. 
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Figure  7.  Calculated  partial  cross  sections  for  the 
excitation  of  the  first  vibrational  level  of  CO  by  collision 
with  Ar  for  a  soft  core  potential  at  a  collision  energy  of 
1.5  eV. 

A,  coupling  retained  between  the  first  two  vibrational  levels; 

B,  coupling  retained  between  the  first  three  vibrational  levels; 
‘  C ,  coupling  retained  between  the  first  four  vibrational  levels. 
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Table  5 


Total  cross  sections  for  the  excitation  of  CO  to  the 
first  vibrational  level  by  Ar  impact 


No.  of  states  coupled 

2 

2 

3 

4 

Potential 

LJ 

SC 

SC 

SC 

* 

Collision  energy 
eV 

Excitation  cross 

2 

section  ) 

** 

0.266 

0.00 

0.00 

0.00 

0.00 

0.600 

0.20~7 

0.13"7 

— 

- —  — 

„  -4 

„  -5 

1.000 

0.35 

0.21 

— - 

— 

1.500 

0.20"2 

0.35"4 

0.54“5 

0.43"5 

2.000 

0 . 14” 1 

0.25“3 

— _ 

—  _  —  _ 

2.500 

0.42"1 

0 , 80-3 

—  —  —  — 

_ 

Cross  sections  for 

Table  6 

the  vibrational  excitation  of 

CO  by  He 

Potential 

LJ 

SC 

SC 

SC 

* 

Collision  energy 
eV 

Excitation  cross 
Q01  Q01 

2 

sections  (na  ) 

°02  °12 

0 . 266 

0.0** 

0.00 

~  _ 

—  — —  — 

0.600 

1.6-2 

1.00"3 

0.77"9 

6.72"6 

1.000 

2.4"1 

1.28"2 

2.60”6 

4.40"3 

1.500 

9.7"1 

1.13"1 

1.66“4 

2.82"2 

2.000 

1.9 

4.06"1 

1.89"3 

8.68"2 

2.500 

2.6 

8.72"1 

8 . 58-3 

1.64"1 

Collision  energy  relative  to  the  ground  state. 

Superscript  indicates  power  of  10  by  which  number 
must  be  multiplied. 
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8.  Discussion 


The  decrease  in  the  QQ1  cross  section  shown  in  tables 
5  and  6  resulting  from  the  use  of  the  soft  core  potential  in 
place  of  the  Lennard-Jones  function  leads  in  turn  to  correspond¬ 
ingly  marked  increases  in  the  calculated  vibrational  relaxation 
times  for  CO  in  the  inert  gases. 

Figure  8  shows  the  theoretical  relaxation  time  for  CO 
in  He  obtained  by  repeating  the  calculation  discussed  in  section 

5  on  the  basis  of  the  Qnl (SC)  cross  sections  given  in  table  6, 

U1  (5) 

in  comparison  with  the  experimental  results  given  by  Millikan 


It  can  be  seen  that  the  agreement  with  experiment  is  now 
much  improved.  The  theoretical  results  lie  within  a  factor  of 
three  of  experiment,  except  at  temperatures  below  300°K  where 
there  is  a  sharp  decrease  in  the  temperature  dependence  of  the 
calculated  relaxation  time. 

In  view  of  the  additional  approximations  that  were  required 
to  establish  the  soft  core  potential  for  the  CO-He  system  this 
result  is  fully  satisfactory  and  would  seem  to  justify  the 
treatment  of  CO-inert  gas  interactions  proposed  in  this  report. 

However,  a  similar  result  is  not  obtained  for  the  CO-Ar 
mixture.  Figure  9  shows  the  relaxation  times  obtained  by  repeating 
the  calculation  discussed  in  section  5  on  the  basis  of  the  two 
state  CO-Ar  QQ1(SC)  cross  sections  given  in  table  5,  in  comparison 
with  the  experimental,  result.  It  can  be  seen  that  the  theoretical 
relaxation  time  is  now  about  a  factor  of  50  greater  than  experiment, 
whereas  the  previous  (LJ)  values  were  about  a  factor  of  4  less. 

It  is  also  apparent  from  table  5  that  coupling  in  additional 
vibrational  states  to  the  calculation  does  not  improve  the  agree¬ 
ment.  Unlike  the  Lennard-Jones  potential,  or  the  CO-He  calculation, 
where  coupling  in  higher  states  in  general  leads  to  an  increase 
in  the  excitation  cross  sections,  for  the  soft  core  potential  a 
decrease  is  produced  at  1.5  eV  which  yields  even  smaller  vibrational 
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relaxation  times.  Figure  7  and  table  5  show  the  result  of 
coupling  three  or  four  states.  Coupling  four  states  produces 
a  further  small  reduction  in  the  cross  sections,  but  these 
results  are  little  different  from  those  for  the  three  state 
calculation. 

This  disagreement  with  experiment  is  about  an  order  of 
magnitude  wor-.e  than  the  poorest  results  obtained  from  the 
previous  applications  of  the  computer  code.  Evidently  the 
discrepancy  is  significant  and  is  all  the  more  surprising  since 
the  CO-Ar  interaction  has  no  peculiar  features  but  should  be  well 
represented  by  the  spherically  symmetric  approximation  used  here. 

As  can  be  seen  the  cross  sections  for  the  excitation  of  CO 
by  Ar  are  very  small  and  the  possibility  that  the  poor  agreement 
with  experiment  is  a  consequence  of  the  numerical  rather  than 
the  physical  treatment  of  the  problem,  possibly  caused  by  the 
weak  coupling  between  the  scattering  equations,  has  been  investigated. 
The  arbitrary  parameters  defining  such  characteristics  as  the 
relative  scaling  between  the  coupled  functions,  the  position  of 
the  starting  point  for  the  step  integration  procedure  and  step 
size  in  the  solution  of  the  scattering  equations  have  been  varied 
over  wide  ranges  without  altering  the  results  more  than  a  few 
percent. 

The  small  cross  sections  are  associated  with  very  extensive 
cancellation  occurring  throughout  the  calculation.  In  fact 
double  precision  on  the  IBM  360  is  insufficient  for  cross  sections 
this  small  for  other  than  nearest  neighbor  transitions  to  be 
reliably  calculated.  To  determine  whether  this  loss  of  precision 
was  introducing  serious  errors  into  the  Qq1  calculation  the  code 
was  run  in  double  precision  on  a  CDC  6600  for  a  collision  energy 
of  1.5  eV.  Although  the  multiple  jump  cross  sections  were  now 
accurately  obtained  it  was  found  that  again  was  virtually 

unaltered. 

Finally,  an  essentially  similar  study  of  ^-Ar  mixtures 
was  carried  out  by  Kindt  and  Marriott  ,  again  on  the 
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This 


basis  of  a  potential  measured  by  Jordan  et  al . 
calculation  simply  confirmed  the  CO-Ar  results.  The  N2  cross 
sections  were  somewhat  larger  and  numerically  better  defined 
than  the  corresponding  values  for  CO  but  the  relaxation  times 
were  nevertheless  within  a  factor  of  two  of  those  for  CO*  as 
should  be  the  case. 

It  can  only  be  concluded  that  the  numerical  results  are 
quite  definite  and  stable  for  the  form  assumed  for  the  soft 
core  scattering  potential  and  that  the  reason  for  disagreement 
with  experiment  must  be  looked  for  elsewhere. 

Because  of  the  low  efficiency  of  Ar  in  vibrationally 
de-exciting  CO  the  experimental  neasurements  will  be  very 
sensitive  to  trace  contamination  by  molecules  such  as  1^0  or 

H0  which  are  very  effective  de-exciters  of  CO.  Using  the 

z  (5) 

relaxation  data  for  CO  in  H2  given  by  Millikan  the  degree 

of  contamination  that  would  be  required  to  account  for  the  dis¬ 
agreement  between  theory  and  experiment  can  be  estimated.  Where 
f  is  the  fractional  contamination  by  H~  the  relaxation  time  t  . 
for  CO  in  the  mixture  will  be  given  by  the  expression 


-1 

T  . 

mix 


(1  -  f) 


-1 

rCO-Ar 


+  f 


-1 

rCO-H 


2* 


At  2000°K,  for  example,  the  computed 
the  measured  relaxation  time  in  Ar, 
and  the  measured 


_ 

1.2 

X 

10  secs 

co-Ar 

T  .  — 

3.5 

X 

-4 

10  secs 

mix 

Tco-H2  = 

7 

X 

_7 

10  secs 

It  follows  that  the  degree  of  contamination  by  H2  that  would 
account  for  the  results  obtained  is  f  «  .002,  or  2000  ppm, 
which  is  two  orders  of  magnitude  larger  than  the  minimum  purity 
level  of  "high  purity  grade  Ar"  supplied  by  the  Matheson  Company. 
Thus  the  lack  of  agreement  also  cannot  be  attributed  to  experimental 
conditions  in  any  simple  way. 

At  this  stage  the  question  remains  unresolved.  More  work 
needs  to  be  done  if  cross  sections  for  heavy  body  collisions  are 
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to  be  calculated  with  any  confidence  for  regimes  where  no 
experimental  work  has  been  performed. 

In  particular  the  possibility  that  collisions  involving 
particles  as  heavy  as  Ar  are  far  more  sensitive  to  the  form  of 
the  scattering  potential  than  is  the  case  for  lighter  atoms 
such  as  He  requires  investigation. 
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Figure  8.  Vibrational  relaxation  time  for  CO  in  He  at  one 

/  c  \  1 

atmosphere  pressure.  Experimental  results  of  Millikan'  . 
Theory  (LJ) ;  calculated  in  this  paper  using  a  Lennard- Jones 
potential  with  coupling  between  up  to  the  first  4  vibrational 
states.  Theory  (SC) ;  calculated  in  this  paper  using  a  soft  . 
core  potential  with  coupling  between  up:  to  6  vibratioral 


states . 
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Ill .  Binary  Gas  Code  and  Molecular  Laser  Modeling 


In  modeling  the  transient  and  steady  state  operation  of 
molecular  laser  systems,  an  understanding  of  vibrational  re¬ 
laxation  processes  in  addition  to  the  effect  of  electrons  and 
inert  diluents  must  be  obtained.  This  Section  of  the  report 
will  describe  the  details  of  the  model  code  development  and  its 
application  to  gas  relaxation  systems. 

The  first  part  of  this  Section  describes  the  relaxation  of 
a  binary  system  of  diatomic  molecules  due  to  both  vibration- 
translation  and  vibration-vibration  processes.  As  an  appli¬ 
cation  of  this  part  of  the  laser  model  code,  calculations  are 
presented  on  the  system  H2  +  D2  2HD.  Following  this,  the 
electron-vibration  rate  coefficients  are  generated  together 
with  spontaneous  emission  coefficients  for  levels  up  to  v  =  29 
in  CO  (both  fundamental  and  overtone) .  These  latter  coefficients 
are  needed  to  compare  the  laser  model  calculations  directly 
to  side-light  measurements  on  laser  systems.  Lastly,  the 
transient  response  of  the  N2-CO-He-e  system  is  shown  as  a 
preliminary  calculation  on  an  actual  molecular  discharge  laser 
system.  In  this  calculation,  30  levels  of  both  N2  and  CO  are 
included  in  the  code. 


III.l  Vibrational  Relaxation  of  a  Binary  Gas  Mixture 

It  is  well  known  that  molecular  species  are  characterized 
by  widely  disparate  vibrational  relaxation  times  in  pure  gas 
systems.^-  The  large  variation  may  be  illustrated  by  the  di¬ 
atomic  molecules  l0,  C1-,  0~,  and  N„  which  have  relaxation 

^  ^  —  P5  £  7  _  5  _  i 

times  at  1000°K  of  5  x  10  ,  5  x  10  ,  7  x  10  ,  and  10  atm-sec., 

2 

respectively.  These  differences  increase  at  lower  temperatures 
while  at  high  temperatures  all  gases  approach  relax  with  unit 
collision  efficienty. 
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In  treating  the  vibrational  relaxation  of  gas  mixtures  in 
circumstances  where  the  values  for  the  pure  gases  are  widely 
separated,  it  is  no  longer  a  straight  forward  problem  to 
correctly  designate  the  relaxation  time.  The  problem  is, 
however,  of  wide  spread  applicability  since  rarely  is  one 
interested  in  knowing  the  behavior  of  pure  systems  (i.e., 
those  systems  composed  on  only  one  species) .  This  is  particu¬ 
larly  true  of  molecular  laser  systems  where  many  species,  both 
atomic  and  molecular  may  be  present.  The  application  of  vi¬ 
brational  relaxation  theories  to  mixtures  is  also  fundamental 
to  many  experimental  techniques  for  determining  vibrational 
temperatures  through  tracer  methods  and  in  determining  the 
effect  of  impurities  on  relaxation  behavior. 

In  general,  there  are  two  important  processes  by  which 
vibrational  energy  is  transferred  in  mixtures.  One  is  the 
process  or  processes  by  which  vibrational  energy  is  degraded 
into  other  forms  of  energy  such  as  translation  and  rotational 
energies.  The  second  process  is  vibration-vibration  exchange 
which  tends  to  couple  the  vibrational  energy  of  the  various 
gases  in  a  mixture.  It  is  the  interrelationship  of  these 
two  processes  which  determines  the  relaxation  of  a  gas  mixture. 

4 

In  earlier  work,  we  presented  an  analysis  and  discussion 
on  the  relative  relaxation  behavior  of  harmonic  and  anharmonic 
oscillators  due  to  vibration-vibration  exchange  processes. 

The  results  of  that  study  can  be  summarized  for  the  case  of 
anharmonic  oscillators  as  follows:  (1)  The  steady  state  vi¬ 
brational  distribution  due  to  vibration  exchange  processes  is 
non-Boltzmann,  in  which  the  upper  levels  are  overpopulated, 

(2)  the  time  to  reach  a  steady  state  distribution  is  an  in¬ 
creasing  function  of  vibrational  level, ^  (3)  total  vibrational 
energy  is  not  conserved  on  the  steady  state  time  scale,  and 
(4)  under  conditions  of  high  vibrational  "temperatures"  and 
low  kinetic  temperature,  population  inversion  in  the  upper 
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vibrational  levels  is  predicted.  These  conclusions  are  markedly 

7  8 

different  from  harmonic  behavior.  ' 

III. 2  Dynamic  Behavior  of  Gas  Mixtures 

4 

In  this  section  we  extend  the  previous  analysis  to  the 
dynamic  behavior  of  relaxing  mixtures  of  diatomic  gases,  and 
herein  discuss  the  various  relaxation  processes  involved  and 
their  respective  relaxation  times  based  on  anharmonic  models. 

4  9 

For  a  pure  gas  it  has  been  previously  established  ' 
that  a  two  time  scale  separation  in  the  relaxation  of  initially 
non-equilibrium  distributions  is  a  valid  assumption  for  the 
lower  levels  of  diatomics  such  as  N2  or  CO  at  low  translational 
temperatures.  A  short  time  scale  is  associated  with  vibration 
exchange  processes  and  a  longer  time  scale  describes  the  re¬ 
laxation  to  equilibrium  through  vibration-translation  degra¬ 
dation  processes.  In  mixtures  of  gases  we  have  a  similar 

4 

situation.  Rapid  vibration  exchange  processes  in  each 
species  will  redistribute  vibrational  energy  into  a  steady 
state  distribution  which  will  slowly  decay  to  equilibrium. 

In  addition  to  these  processes,  however,  the  relaxing  distri¬ 
butions  will  be  affected  by  vibrational  coupling  between  the 
species  through  vibration  exchange  reactions. 

We  have  previously  shown  the  strong  dependence  of  vibra¬ 
tional  exchange  rates  on  energy  defect‘d  (i.e.,  the  energy 
transferred  into  or  out  of  relative  translational  energy  during 
an  exchange)  and  consequently,  the  coupling  between  two  species 
in  a  mixture  is  expected  to  be  strongly  dependent  on  the  re¬ 
spective  vibrational  spacing.  For  mixtures  such  as  N2  and  CO 
where  the  vibrational  spacing  of  the  gases  is  similar  (i.e., 

E,  =  188  cm-1  for  the  v.n  -  vnl  levels  and  E  =  14  cm-1  for  the 
d  01  01  d 

vcn  -  v-t  levels  of  N0  and  CO,  respectively)  ,  the  exchange 

b  /  ui.  c. 

rates  within  each  species  will  not  display  an  energy  defect 
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while  rates  between  species  will  be  characterized  by  a  rather 
large  defect  (i.e.,  100  cm  ^)  ,  thus,  in  the  harmonic  model,  a 
slight  time  separation  in  the  pure  gas  vibrational  exchange 
processes  and  the  coupling  between  species  is  expected. 

For  a  mixture  of  gases  such  as  N2  and  NO,  the  energy 
defect  between  species  is  considerably  larger  than  that  ex¬ 
hibited  in  either  of  the  pure  gases  and  consequently,  three 
distinct  time  scales  should  be  evident.  The  first  time  scale 

will  feature  rapid  vibrational  rearrangement  in  each  molecule 

'  4 

leading  to  a  steady  state  distribution  as  previously  predicted. 
A  second  and  longer  time  scale  will  describe  the  coupling  of  the 
two  vibrational  distributions  to  form  a  second  steady  state 
distribution  which  will  then  proceed  toward  equilibrium  through 
vibration-translation  processes  on  a  third  time  scale. 

As  the  energy  defect  between  two  species  in  a  mixture  is 
increased  by  suitably  changing  one  of  the  species,  eventually 
a  merging  of  the  time  scales  for  mixture  coupling  and  vibra¬ 
tion-translation  processes  for  one  of  the  species  will  occur. 
This  situation  is  evident  for  N2  and  02  mixtures  at  low  tem¬ 
peratures  (i.e.,  <  3000°K)  as  demonstrated  by  the  measurements 

3 

of  Millikan  and  White. 

In  order  to  illustrate  these  various  time  scales  for  the 
case  of  mixture  relaxation,  the  master  equation  for  the  time 
evolution  of  a  binary  mixture  has  been  formulated  and  numer¬ 
ically  integrated. 


Ill 
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The  Master  Equation 


The  master  equation  for  a  binary  mixture 


may 


be  written: 


4,7 


dn. 


dt 
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£  £ 
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SR 

(Q  n  N  -  Q 
'  sr  s  s 


RSn  N  ) 
rs  r  R 


0) 


in  which  n.  is  the  population  of  the  i-th  vibrational  level; 
n  is  the  total  number  density  of  species  A;  and  N,  the  cor¬ 
responding  densities  of  species  B;  Pgr  is  the  rate  constant  for 
the  process: 


n„  +  n  — ►  n„  +  n 
s  r 


(2) 


in  which  vibrational  energy  corresponding  to  -  Eg,  the 
difference  between  the  energies  of  the  r-th  and  s-th  states 
in  species  A,  is  transferred  to  or  from  translational  energy; 

1  m 

P  is  the  rate  constant  for  the  process: 

sr  * 


n 


+  n. 


n  +  n 
r  m 


(3) 


in  which  vibrational  energy  is  exchanged  within  species  A, 

the  slight  energy  defect  being  compensated  by  translational 

SR 

energy;  and  Q  and  Q  are  the  analogous  rate  constants 

SL  SL 

for  the  vibration-translation  and  vibration-vibration  energy 
transfer  processes  involving  the  cross  terms  between  species 
A  and  B: 


and 


n„  +  N  — ►  n_.  +  N 
s  r 


(4) 
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(5) 


ns  +  Ns  — ’  nr  +  nr  ' 

4 

In  our  earlier  work,  we  concluded  that  in  the  pure  gas 
case,  only  single  quantum  processes  need  be  considered.  How¬ 
ever,  for  mixtures  of  gases  with  greatly  differing  vibrational 
spacing  there  is  a  possibility  that  resonant  vibration- vibra¬ 
tion  exchange  might  involve  a  double  quantum  transition  in 
one  of  the  species.  in  general,  for  the  lower  vibrational 
levels,  these  processes  will  be  slow  in  comparison  to  single 
quantum  processes  due  to  the  form  of  the  rate  expressions^® 
and  we  will  restrict  ourselves  in  this  report  to  systems 
wherein  only  single  quantum  processes  are  important.  Evi¬ 
dently  this  is  not  a  major  limitation  since  most  of  the 
diatomic  systems  of  interest  fall  within  this  classification. 

This  restriction,  coupled  with  the  use  of  the  detailed 
balance  conditions. 


-(Er  +  1 

-  Er)/kT 

Pr,  r  + 

1  "  Pr 

+ 

1,  r  e 

4  1 

l  =  Qr 

“(Er  +  1 

-  /kT  i 

^r,  r  + 

+ 

1,  re 

ml 

Prs  = 

Plme 
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(Es  +  Ej  -  Er 

-  E  )/kT 
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and 


RS  ^  SR 
Q  =  Q  e 


■(Es  +  F  -  Er  -  F  )/kT, 


‘rs 


sr 


where  Fj  represents  the  vibrationel  energy  of  the  I-th  level 
of  species  B  above  the  zeroth  level,  permits  the  master  equa¬ 
tion  to  be  written  as: 
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(7) 


where  &r  is  a  production  or  destruction  source  of  nr.  In  later 
sections  of  this  report,  this  term  will  include  the  excitation 
of  vibrational  levels  by  hot  electrons  and  the  spontaneous  emis¬ 
sion  associated  with  infrared  active  molecules,  e.g.  CO. 

Equation  (7)  formally  describes  the  rate  of  change  of  nr 
with  time  in  the  absence  of  dissociation.  The  effect  of  disso¬ 
ciation  on  vibrational  distributions  is  currently  under  study. 
Three  distinct  types  of  terms  are  present.  The  most  familiar 
are  those  for  the  transfer  of  vibrational  to  translational 
energy  and  the  reverse  (i.e.,  P  . n  n) .  These  processes  will 

”  S  p  S  •  X  s 

be  designated  as  VT.  The  vibration- vibration  exchange  terms  in 
a  single  species  (i.e.,  ^^+1* r  nr+l  ns-l^  be  designated  by 

WA  A  or  WB_B,  while  the  exchange  terms  between  species 
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(i-e"  “hI^Vr-I1  Wil1  be  n°ted  as  WA-B- 

The  set  of  equations  (7)  with  the  rate  constants  outlined 

in  the  next  Section  have  been  numerically  integrated  by  the 

14 

Runge-Kutta-Merson  technique.  This  numerical  technique  has 

been  used  extensively  for  the  integration  of  large  sets  of 

coupled  non-linear  ordinary  differential  equations  such  as  are 

found  in  modelling  coupled  chemical  reactions.  A  particular 

14 

feature  of  the  code  as  taken  from  the  work  of  Keneshea  is 
the  use  of  a  version  of  the  steady  state  (S-S)  approximation 
for  those  equations  which  satisfy  an  appropriate  criterion. 

In  a  coupled  set  of  rate  equations  such  as  occur  for  vibrational 
relaxation,  the  equation  can  be  written 


=  IX  -IX  •  nr  <8> 

where  the  F  and  ^Rr  terms  may  be  complicated  functions 
of  all  the  other  concentrations  and  rate  coefficients.  Under 
conditions  where  the  magnitudes  of  both  ^Fr  and  the  term 
Rr  •  nr  are  large  but  their  difference  very  small,  the 
time  progress  in  the  normal  Runge-Kutta  integration  is  very 
slow,  i.e.  the  set  of  equations  become  "stiff."  In  the  Keneshea 
version  of  the  code  as  used  here,  each  equation  is  checked  to 
determine  if  it  satisfies  the  S-S  criterion,  meaning  the  mag¬ 
nitude  of  the  ^Fr  term  is  large  but  the  derivative  is  small. 
When  this  criterion  is  passed  the  equation  is  iterated  over 
the  At  time  step  using  an  algebraic  algorithm  as  shown  below. 


For  conditions  under  which  ^F^  and  ^Rr  are  both  large, 
rather  than  use  the  canonical  S-S  assumption 


dn 

_ _ i 

dt 


=  0 


(9) 


which  leads  to 


III 


8 


(10) 


n 


r 


the  more  realistic  assumption  is  made  that  the  and 

/  iRr.  terms  are  effectively  constant  over  the  At  step  leading  • 
to  a  first  order  ordinary  differential  equation.  The  solution  is 


nr (t  +  At) 


(11) 


This  solution  is  used  for  all  equations  satisfying  the  steady 
state  criterion  and  is  iterated  over  a  At  step  together  with 
the  other  equations  which  continue  to  be  integrated  by  the 
Runge-Kutta-Merson  technique.  When  an  equation  ceases  to 
satisfy  S-S,  it  is  automatically  returned  to  the  R-K-M  sub¬ 
routine  for  integration. 


By  the  use  of  this  S-S  algorithm,  the  time  step  in  the 
integration  can  be  kept  large  while  not  sacrificing  either 
stability  or  accuracy  in  the  coupled  set  of  equations. 


The  vibration-vibration  rate  coefficients  used  in  the 

code  are  adapted  from  a  model  presented  by  Rapp  and  Englander- 
(12) 

Golden  .  The  derivation  of  the  rate  coefficient  is  given 
in  Appendix  I  while  the  actual  expression  used  is 


Prn  =  1.508-12,7/6L10/3a2E//3T5/6U2  U2 
sm  d  rn  sm 


where 


•  exp  -1.32 (LEd) 2/3  $  1/3 


(12) 


ix  =  reduced  mass  in  molar  units 

0 

L  =  exponential  potential  parameter  (A) 

0 

n  =  hard  sphere  diameter  (A) 

T  =  kinetic  temperature  (°K) 


Ed  =  energy  defect  between  vibrational  quanta  being 
exchanged  (cm_l) 
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and 

U  =  matrix  element  for  the  r  -  n  vibrational  transition 
rn 

The  matrix  elements  used  are  modified  Harmonic  Oscillator 
(13) 

matrix  elements'  in  which  the  actual  energy  spacing  between 
vibrational  units  is  used  instead  of  the  constant  spacing 

assumed  in  the  H-D  model.  Very  close  agreement  has  been  fqund 

! 

for  these  matrix  elements  by  comparison  to  Morse  oscillator, 
results  as  shown  in  Table  I. 

Vibration-translation  rate  coefficients  have  been  taken 
from  measured  relaxation  times  using  a  harmonic  oscillator 
model  The  expression  used  is 

P  =  *1 

rn  pr 

where  pT  is  the  measured  relaxation  time.  Scaling  with  vi-  * 
brational  level  is  given  by  the  square  of  the  vibrational 
matrix  elements.  , 

1 1 1 . 3  Application  of  the  Binary  Code  to  the  Reaction 

As  a  test  for  the  model  code  and  to  interpret  the  results 
of  recent  experiments  on  the  reaction  , 

H2  +  D2  t  2HD  ,  (14) 

the  vibrational  relaxation  code  has  been  applied  to  this 
system.  Previous  studies ^ have  noted  the  following  facts 
regarding  this  reaction  : 

(a)  in  the  phenomenological  rate  expression,  the  ex¬ 
ponent  of  the  concentrations  of  H2  and  D2  are  not 
the  same, 

(b)  the  rate  appears  quadratic  in  time, 

(c)  vibrational  energy  in  either  H2  or  D2  appears 
necessary  to  initiate  the  reaction. 


1  -  e 


-hcw/kT\- 


(13) 
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I  I 


Table  i 

BOUND-BOUND  MATRIX  ELEMENTS 
FOR  N.  ' 


Element 


Harmonic 


Morse 


0-1 

1-2 

2- 3 

3- 4 

4- 5 

5- 6 

6- 7 
,7-8 

8- 9 

9- 10 
10-11 
11-12 

12- 13 

13- 14 

14- 15 
'15-16 

16- 17 

17- 18 

18- 19 

19- 20 

20- 21 
21-22 

22- 23 

23- 24 

24- 25 

25- 26 

26- 27 

27- 28 

28- 29 

29- 30 

30- 31 

31- 32 

32- 33 

33- 34 


.08058 
.11505 
.14227, 
.161589 
.  19729 
.20721 
.22607 
.24413 
.26160 
.27862 
.29529 
.31169 
.32790 
.34397 
.35996 
.37589 
.39182 
.40777 
.42376 
.43,984 
.45603 
.47234 
.48881 
.50545 
.52229 
.53935 
.55665 
.57422. 
.59208 
.61023 
.62872 
.64756 
.66680 
.68642 


.08062 
.11366, 
.13876 
.15973 
.17800 
.19441 
.20927 
.22305 
.23580 
.24771 
,25900 
.26959 
.27970 
' .28927 
.29854 
.30716 
.31575 
.32363 
.33154 
,33889 
.34597 
.35298 
.35964 
.36623 
.37228 
,  .37855 
.38445 
,.39014 
.39541 
.40099 
.40592 
.41108 
.41618 
.42109 
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Theoretical  studies  have  not  been  successful  in  under¬ 
standing  the  small  activation  energy,  i.e.  42  kcal/mole, 

which  is  considerably  less  than  the  dissociation  energy  of 

(17) 

either  Hj  or  Dj,  although  studies  of  this  4  center 

system  cannot  be  fully  treated. 

The  particular  experiments  modeled  by  detailed  calcula¬ 
tions  involve  the  Raman  excitation  of  the  v  =  1  level  of  Hj 
in  the  presence  of  D2  and  the  subsequent  formation  of  HD. 

The  reaction  vessel  initially  has  5  atm  each  of  Hj  and  D2 
at  room  temperature  and  the  raman  excitation  is  postulated 
to  give  roughly  1%  H2(v  =  l)  per  pulse.  Following  many 

pulses,  the  reaction  mixture  is  analyzed  for  HD.  The  exper- 

+16 

iments  suggest  approximately  5  x 10  HD  molecules/cc  produced 
per  pulse. 

In  performing  this  model  calculation,  in  addition  to 
the  W  and  VT  processes  which  are  operative  in  and  between 
H2  and  D2,  the  following  reaction  terms  were  included  in 
the  code: 


icT12 

for 

v  a  3 

H-(v) 

+  D-  t- 

2  HD 

k(v)  = 

(15) 

4 

4 

0 

for 

v  <  3 

|10-12 

for 

v  2  4 

D-(v) 

+ 

2  HD 

k  (v)  = 

(16) 

4 

4 

1 0 

for 

v  <  4 

The  rate  coefficients  are  assumed  to  be  zero  until  a  vibra¬ 
tional  level  is  reached  approximately  consistent  with  the 

measured  activation  energy.  The  VT  rate  coefficients  used, 

(18) 

i.e.  the  relaxation  times  ,are  taken  from  Kiefer  and  Lutz' 
while  the  W  rate  coefficients  are  calculated  from  the  model 
in  the  previous  section. 

In  Figure  1  is  shown  the  transient  behavior  of  represen¬ 
tative  vibrational  levels  and  in  Figure  2,  the  HD  produced  via 
the  preceeding  reactions  for  the  case  in  which  1%  H2(v  =  l)  is 
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TIME  (SEC) 

Figure  1A:  Relaxation  of  Vibrational  Energy  in  H_-  D_  Mixtures. 


Figure  2  -  HD  Production  in  H_  +  D_  +-  2HD 


excited.  The  results  show  production  rates  of  HD  consistent 
with  experimental  observation  but  the  most  important  point 
is  that  even  though  vibrational  levels  where  initially 
excited,  the  majority  of  the  HD  formed  is  through  vibrationally 
excited  D2.  This  is  a  direct  result  of  the  W  processes  which 
tend  to  preferentially  pump  D2  levels  from  H2  due  to  the  smaller 
vibrational  spacing.  The  time  scales  for  this  exchange  are 
displayed  in  Figure  1.  It  is  important  to  note  that  if  the 
initial  conditions  are  1 %  D~(v  =  l),  then  the  HD  produced  is 

^1  f. 

again  of  the  order  of  5  x 10  molecules/cc,  but  for  this  case 

there  is  essentially  no  contribution  due  to  the  vibrationally 

excited  H2  reaction,  as  expected  from  the  larger  vibrational 

spacing.  The  details  of  this  calculation  are  currently  being 
(19) 

prepared. 


I I I . 4  Excitation  of  and  CO  Vibrational  Levels  by  Electrons 

The  excitation  of  vibrational  levels  of  N2  and  CO  by 
electrons  is  an  important  process  in  molecular  discharges, 
molecular  lasers  and  in  the  upper  atmosphere.  In  studying 
the  transient  and  steady  state  behavior  of  molecular  systems 
in  which  hot  electrons  are  present,  the  rate  coefficients 
characterizing  the  coupling  between  the  electrons  and  vibra¬ 
tional  levels  are  required.  In  this  section,  rate  coefficients 
are  calculated  from  cross  section  data  of  Schulz and 
Chen'  .  Of  particular  interest  is  the  temperature  depen¬ 
dence  of  the  rate  coefficients  at  elevated  temperatures 
(appropriate  to  molecular  laser  systems)  and  the  scaling  of 
the  rate  coefficients  with  increasing  vibrational  level. 


Ill 
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Analytical  Development 


The  rate  coefficient,  k(T)  can  be  defined  in  terms  of 
the  cross  section,  Q(v) ,  as 


o& 

k(T)  =|  Q(v)  v  f  (v)  dv  (17) 

“  «0 


where  f(v)  is  the  Maxwell-Boltzman  distribution  of  velocities. 
By  transforming  f(v)dv  to  a  distribution  of  speeds,  equation 
(17)  can  be  written 


k(T) 


4tt 


v3 


exp 


-mva 

2kT 


dv. 


(18) 


Equation  (18)  can  be  easily  changed  to  the  energy  represen 
tation,  yielding  the  following  rate  coefficient: 


k(T) 


8rr 


3/2 

1  1/2  r 

( ,  i 

( 

— J  Q(E)  E  exp 

“E/kT 

\  2rrkT 

)  \ 

»  j  0 

2 

where  E  =  ^m v  ,  m  being  the  reduced  mass  of  the  collision 
(essentially  the  electron  mass)  . 


(19) 


The  available  data  on  electron  excitation  cross  sec¬ 
tions  for  N2  and  CO  can  be  fit  over  a  limited  electron 

energy  range  E^  -  E2 ,  where  Ej  >  E^  by 


Q(E)  =  AE  +  B 


(20) 


where  A  and  B  are  constants  determined  by  the  end  point 
values  E^ ,  and  E2 . Qj •  Thus 
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Q8  -  Qi 


E2  -  Bx 

EaQi  "  *2^2 

i7rEi 


(21) 


Over  the  energy  range  -  E2  =  AE12,  the  integral  in 
(19)  can  be  integrated,  yielding 

E.  -3/2  -1/, 


k(T)  |  =  8rr(2rrkT)  ~  (m) 

e2 

E, 


+AE2 


and 


E, 


k  (T)  =  ^  MT) 


Z-J 

A  El? 


*~E/kT 

e  /XT  (kT) 


(2kTA+B)  (E+kT) 


(22) 


E, 


E. 


(23) 


Equations  (22)  and  (23)  have  been  programmed  to  yield 
excitation  rate  coefficients  as  a  function  of  the  electron 
temperature  using  input  cross  section  data. 

Results  and  Discussion 


The  rate  coefficients  for  the  excitation  of  vibrational 

levels  of  N2  and  CO  by  electrons  in  the  electron  temperature 

range  1000°K  to  30,000°K  using  the  measured  cross  section 

data  of  Schulz are  shown  in  Figures  3  and  4,  respectively. 

The  complete  numerical  results  for  these  cases  are  available 
(22) 

elsewhere  ,  together  with  the  cross  section  data  used  in 
the  calculation.  In  general,  the  error  limits  on  the  calcu- 
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VIBRATIONAL  EXCITATION  RATE  COEFFICIENTS 
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lated  rate  coefficients  are  about  ±20%. 

The  spacing  of  the  electron  excitation  rate  coefficients 
for  progressively  larger  quantum  changes  reflects  the  negative 
ion  mechanism  responsible  for  the  excitation  process.  In 

f21  23) 

both  N0  and  CO,  the  negative  ions  are  unstable  by  1.89  ev'  '  ' 

1  (21) 

and  1.3  ev  ,  respectively,  corresponding  to  about  the  6th 
vibrational  level  of  each  molecule.  The  spacing  of  the  rate 
coefficients  indicates  a  rather  close  coupling  between  levels 
near  v  =  6  while  quantum  changes  greater  than  6  show  a  con¬ 
siderably  reduced  coupling,  i.e.  a  wider  spacing. 

The  excitation  rate  coefficients  calculated  from  the 

( 2 ) 

theoretical  cross  sections  of  Chen'  '  are  tabulated  together 

(22) 

with  the  corresponding  cross  sections  elsewhere  .  The 
general  agreement  between  the  calculated  and  measured  values 
is  within  the  accuracy  of  the  overall  rate  coefficients, 
i.e.  20%. 

By  observing  the  calculated  rate  coefficients  derived 
from  Chen's  cross  sections,  some  indication  of  the  scaling 
necessary  with  increasing  vibrational  level  can  be  obtained. 
Figures  5  and  t>  show  this  relationship  for  Av  =  1  and  Av  =  5 
processes,  respectively.  No  trend  in  the  calculated  values 
are  noted  with  increasing  vibrational  level  in  the  electron 
temperature  range  between  10,000°K  and  30,000°K,  therefore, 
to  within  the  accuracy  of  these  calculations  all  excitation 
processes  with  the  same  av  change  can  be  assigned  the  same 
value  as  the  v  =  0  excitation  case  at  any  given  electron 
temperature. 

It  is  also  very  important  to  note  that  above  about 
10,000°K  the  electron  excitation  rate  coefficients  are 
effectively  independent  cf  temperature.  This  fact  means  that 
in  calculating  the  vibrational  distributions  in  a  discharge 
environment  in  which  the  electron  temperature  is  greater  than 
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this  value,  no  account  has  to  be  taken  of  the  non- Boltzmann 

(25} 

electron  distributions  calculated  to  be  present'  ‘  .  This 
is  not  the  case  for  excitation  of  electronic  states  since 
these  processes  will  be  very  sensitive  to  the  high  energy 
electron  tail. 


I I I . 5  Spontaneous  Emission  Coefficients  for  CO 


As  part  of  the  modeling  of  molecular  lasers,  the  spon¬ 
taneous  emission  emitting  from  vibrational  transitions  in  the 
infrared  must  be  characterized.  This  serves  two  purposes; 
firstly,  the  emission  from  higher  lying  vibrational  levels 
can  be  a  significant  energy  loss  mechanism  in  distorting  the 
vibrational  distribution  and  secondly,  in  evaluating  the 
model  laser  code  through  side-light  emission  the  code  must 

be  able  to  generate  artificial  spectra.  In  this  section,  we 
present  the  fundamental  and  first  overtone  Einstein  Coefficients 

for  Spontaneous  Emission  calculated  from  the  data  of  Young  and 

„  .  (24) 

Eachus' 


Young  and  Eachus  have  evaluated  4  cubic  dipole  moment 
functions  and  have  concluded  that 

M(r)  =  ±  0.112  +  3.11(r  -  r  )  -  0. 15 (r  -  r  ) 2  -  2. 36 (r  -  r  ) 3 

c  c  c 

(24) 

is  best, based  on  the  relative  rotational  line  intensities  in 
the  overtone  bands  of  CO.  Using  the  integrated  band  inten¬ 
sities  for  the  fundamental  and  overtone  bands  of 


-1  -2 

a,  ,  .  .  =  261  amagat  cm 

fundamental  3 


and 


-1  -2 

a  .  =1.88  amagat  cm 

^overtone 


III  -  23 


the  following  table  of  Einstein  Coefficients  has  been  generated 
and  incorporated  into  the  laser  model  code.  Current  work  is 
concentrated  on  a  subroutine  to  generate  the  entire  artificial 
spectrum  in  the  vibration-rotation  bands  including  an  appro¬ 
priate  slit  function. 


III. 6  Preliminary  Laser  Model  Calculations 

The  several  sections  of  the  numerical  code  as  described 
in  the  previous  pages  has  been  incorporated  into  a  preliminary 
laser  code.  Currently  the  code  is  capable  of  handling  up  to 
40  levels  in  each  of  two  diatomic  molecules,  i.e.  Nj  and  CO, 
while  also  including  excitation  via  hot  electrons  and  VT 
processes  due  to  an  inert  diluent.  The  code  is  particularly 
suited  to  study  the  transient  response  of  the  vibrational 
distributions  of  both  diatomics  to  a  steady  state  electron 
temperature  such  as  would  occur  in  a  direct  discharge  laser. 
Preliminary  calculations  have  been  performed  to  check  out 
the  operation  of  the  code. 

Figure  7  shows  the  time  development  in  the  vibrational 
distribution  of  CO  in  the  presence  of  20,000°K  electrons 
using  a  20  level  model  for  CO  in  the  presence  of  He.  The 
structured  appearance  to  the  early  time  profile  is  due  to 
the  non-uniform  electron-vibrational  coupling  rates  as 
shown  in  Figures  3  and  4.  As  time  progresses,  VV  processes 
become  important  and  smooth  out  the  vibrational  profile. 
Further,  due  to  the  decreasing  vibrational  spacing  with  in¬ 
creasing  quantum  number,  higher  lying  vibrational  levels  are 

pumped  preferentially  at  the  expense  of  lower  ones  as  has 

(4  7  9) 

been  pointed  out  before  '  '  .  At  long  times,  when  the 

lower  vibrational  levels  have  reached  a  steady  state  value, 
the  upper  levels  continue  to  grow  and,  in  fact,  show  popula¬ 
tion  inversion.  This  behavior  of  the  20  level  code  has  been 
traced  to  the  termination  of  the  set  of  vibrational  equa- 
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Table  II 


Relative  Einstein  Coefficients  for  Spontaneous  Emission  for  CO 
upper  vibrational  level  fundamental  overtone 


1 

1.00 

— 

2 

1.93 

1.00 

3 

2.78 

3.00 

4 

3.53 

5.96 

5 

4.21 

9.89 

6 

4.81 

14.74 

7 

5.35 

20.48 

8 

5.74 

27.03 

9 

6.23 

34.36 

10 

6.58 

42.43 

11 

6.88 

51.17 

12 

7.14 

60.54 

13 

7.35 

70.79 

14 

7.60 

81.62 

15 

7.73 

92.99 

16 

7.83 

104.82 

17 

7.90 

117.05 

18 

7.93 

131.84 

19 

7.94 

147.24 

20 

7.83 

160.27 

21 

7.69 

173.80 

22 

7.53 

187.53 

23 

7.35 

201.19 

24 

7.16 

214.94 

25 

7.00 

228.60 

26 

6.74 

242.13 

27 

6.51 

255.45 

28 

6.27 

268.52 

29 

5.82 

281.40 
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tions  at  20.  A  30  level  set  gives  identical  results  for  the 
lower  levels  but  shifts  the  population  to  higher  levels  than 
shown  in  Figure  7.  This  effect  is  under  current  study  as  it 
appears  a  numerical  result  rather  than  one  which  is  physically 
expected. 

The  ciicle  points  in  Figure  7  correspond  to  a  steady 

/  2  g  \ 

state  analysis  performed  by  Rich'  ' .  The  agreement  over  the 
entire  range  is  obtained  in  the  code  with  increased  vibra¬ 
tional  levels. 

The  laser  code  has  also  been  run  with  concentrations  of 
N2  present  but  the  effect  on  the  time  development  of  the  CO 
vibrational  distribution  is  small  for  concentrations  of  N2 
equal  to  or  less  than  those  used  for  CO  due  to  the  decreased 
electron- vibrational  rate  coefficients  for  Nj.  This  particu¬ 
lar  effect  will  be  further  quantified  in  continued  numerical 
studies . 

Lastly,  a  subroutine  is  currently  being  debugged  which 
will  provide  an  artificial  spectrum  of  the  time  developing 
CO  vibrational  distribution  for  direct  comparison  to  time 
resolved  side-light  measurements  on  a  CO-^-He  molecular 
laser.  The  initial  phases  of  this  experimental  program  are 
outlined  in  the  following  section. 
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IV.  Experimental  Facilities 


The  CO  laser  exhibits  a  tremendous  diversity  of  operating 

characteristics.  Laser  action  occurs  either  pulsed ^  or  CW^\ 

at  liquid  N0  temperature^,  dry  ice  temperature^,  or  room 
Z  (5) 

temperature  .  Laser  emission  has  been  observed  on  transitions 
from  v  =  4tov=3,  at  4.9  p,m,  to  v  =  16  to  v  =  15,  at  5.7  y,m. 

In  general,  the  addition  of  various  heavy  gases,  such  as  Hg  or 
Xe,  can  increase  the  efficiency  of  laser  operation. 

In  order  to  monitor  the  energy  transfer  processes  in  a 
mixture  of  CO,  N2,  He,  02  and  other  additives,  an  experimental 
facility  has  been  constructed  as  follows.  An  optical  cavity 
has  been  designed,  consisting  of  two  gold-coated  mirrors,  one 
flat,  and  one  with  a  4  meter  radius  of  curvature.  A  one  mm 
hole  was  drilled  in  the  curved  mirror  for  output  coupling. 

A  one-meter  long,  pyrex  tube  was  constructed  with  four 
viewing  ports  equally  spaced  along  the  upper  surface,  20  cm 
apart.  A  NaCl  window  is  secured  to  each  viewing  port.  The 
entire  discharge  tube  is  enclosed  in  a  plexiglass  cooling 
enclosure.  The  design  of  the  system  is  such  that  an  IRTRAN 
lens  of  f/2  focal  ratio  can  be  used  for  viewing  the  spontaneous 
emission  from  the  discharge. 

At  each  end  of  the  tube  are  end  assemblies  consisting  of 
cylindrical,  hollow  copper  electrodes,  and  Brewster- angled 
NaCl  windows.  A  20kv,  50mA  power  supply  is  used  to  excite  the 
tube,  either  pulsed  or  CW. 

Emission  spectra  are  recorded  using  a  Jarrell-Jvsh  one- 
quarter  meter  infrared  monochromator,  and  either  a  lead  selenide 
or  liquid-N2  cooled  Hg-doped,  germanium  detector,  depending  on 
the  spectral  range  of  observation.  For  CW  observation,  a  PAR 
chopper  and  lock-in  amplifier  are  used  to  provide  phase-sensitive 
detection. 


IV-1 


Synthetic  spectra  are  being  generated  by  computation  to 
guide  the  identification  of  discharge  conditions.  An  improved 
power  supply  has  been  ordered  along  with  a  coherent  radiation 
model  201  power  meter,  to  augment  the  laboratory  capability. 

A  gas-handling  system  capable  of  feeding  up  to  4  gases 
into  the  laser  tube  simulatneously  has  been  constructed.  Flow 
rates  are  monitored  for  each  gas  with  individual  rotometer 
readings.  The  total  gas  pressure  in  the  discharge  tube  is 
continually  monitored. 

In  summary,  a  flexible  experimental  facility  has  been 
designed  and  constructed  for  the  observation  of  both  spon¬ 
taneous  and  stimulated  emission  from  gas  mixtures  relevant 
to  CO  laser  operation,  under  a  variety  of  excitation  condi¬ 
tions.  Survey  spectra  have  been  taken,  and  synthetic  spectra 
are  being  generated  for  identification  of  discharge  conditions. 
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Appendix  1 


VIBRATION  -  VIBRATION  RATE  COEFFICIENTS 


Following  Rapp  (Al)  we  consider  the  following  model. 


We  define 

ilr * B  (xa)  -  wave  function  of  molecule  AB  in  vib.  state  r. 

'I'*0  (x2)  -  wave  function  of  molecule  CD  in  vib.  state  s. 

The  total  wave  function  for  the  system  of  two  interacting  molecules 
can  be  written  as  an  expansion  in  the  time  dependent  wave  functions 
for  each  molecule. 

V  '  “iuu,  t  -iQ.t 

(1)  4  =  L  cri  (+n;B  (X,  H.CD  (x3)  e  e 

r  ■ 


U>r 


It 

h 


and 


where 


We  set  the  initial  conditions  as  the  following: 


(r, A0  (Xj)  =  tlrp  AB  (xx ) 


therefore  Cpq  (-») 


and  i];iCD(x2)  =  i!rqc0  (x2)  as  t  -  -® 

=  1  and  Cr f  (-»)  =  0  for  r  ±  p  and  s  ±  q. 


We  are  now  interested  in  solving  the  Schrodinger  equation  for  the 
coefficient  CBB(+»).  The  Hamiltonian  is  written 

I 

H  =  H* B  (Xj)  +  Hc  0  (xa)  +  V(y,x,,x2) 


where 

HAb(x1)  =  Hamiltonian  for  the  isolated  AB  diatomic, 

Hc 0 (x2 )  =  Hamiltonian  for  the  isolated  CD  diatomic, 

V(y,Xj,xa)  =  Interaction  perturbation. 

The  assumption  is  now  made  that  the  interaction  is  exponentially 
repulsive  in  the  distance  B-C,  and  the  result  obtained  (see  Rapp 
and  Sharp  (Al) )  is 

(3)  V(Y.  Xl(X.)  -E.  exp  j-  £  |y  -{jj^ }  X,  |  . 

We  note  that  this  form  allows  the  interaction  energy  to  be  written 
as 

(4)  V(Y,Xj  ,X2)  =  V(Y)V(X1)V(Xa) 


where 

V(Y)  is  determined  from  the  classical  trajectory  of  the 
relative  motion  and 

Y  =  y-yT  ,  x1=  Xi-x1T  f  x3  =  x2-x2T 


where  the  t  subscripted  variables  represent  turning  point  variables. 


When  the  classical  equation  of  motion  is  used  (A2) ,  one  gets 


where 


V  (Y) 


Eo  =  i^vo  and 


= 


MA  b  Mco 
Mab  +  Mcd  * 


Inserting  the  potential  into  the  Schrodinger  equation,  multiplying 
by  i|r*B(x,),  *«0(xa)  and  integrating  over  xx  and  x2  we  find 


(6) 


Cb«  =  ^i)~1  Eo  sech2 


«  (t)ei(u,»‘u,')t  el(n*  “  fl«)t 

'  '  r  '  i 


*U  U 

“n  •  ■ 


where 


(7) 


00 

UPB  =  J  ilfj0  (XjJdXj  ,  etc 


We  now  use  the  2-state  approximation  and  arrive  at  the  coupled 
equations  of  Rosen-Zener  (Al,  A3) . 

.  -2niAvt 

Cx  =  ife  Cs 

•  +  2-rriAvt 

Ca=  ife  Cj 


which  has  as  a  solution 
(8)  |C1  (“)  I  2  =  sin3  A  sech3  (ttAv  t  ) 


where 


A  =  |  f (t) dt 
•00 


and 


E»  ?  f 

f(t)  *  ~  u-u”  J  seoh  v 

*  00  \ 


2TTV0t 


dt 


8L 


3 


This  can  be  further  reduced  and  integrated  to  give 
4LE„  * 

f  (t)  =  JsechOOdx 

00 

4LE0 


and  therefore 

4n|jLv0 

A - h — :  uPBuiB  . 


(9) 


(10) 


Inserting  this  into  (8) ,  we  obtain 

f  4ttLv0 


PJJ  (v)  =  sin3j- 


UP  BU, ,  (sech3 


where 

H  =  reduced  collision  mass 
L  =  exponential  parameter 
Ea  =  energy  defect 
v0=  relative  velocity 
h  =  Planck's  constant 


4 


and  Pjjj  (v)  =  velocity  dependent  probability  for  vibrational 
exchange  from  r-*n  and  s-*m. 


We  now  consider  two  limits  on  (10) ,  the  case  of  large  energy 
defect  and  small  energy  defect.  First  we  note  for  essentially 
all  cases  of  interest 


(ID 


sin* 


(4tt|j,Lv0 


16tt3 


|j,8L8v0  3 


h3 


U; 


•  ■ 


And,  for  the  case  of  large  energy  defect 


(12) 


sech3 


2LE„ 


v. 


S  4  exp! 


4LEd  > 


) 


4 


(13) 


and  for  the  case  of  small  energy  defect 

/ 2LE4  \  4LaEds 

sechs(  — - )  ^  1  - 


\  vo 


v. 


The  rate  coefficient  for  vibration-vibration  exchange  processes 
can  be  obtained  from  (10)  by  integrating  over  the  collision  fre¬ 
quency  and  the  one-dimensional  velocity  distribution.  Thus, 


(14) 


00 

r;:(t)  =  J  p.r;(v0)  nv0)f(v0)dvc 


where  r(v0)  =  hard  sphere  collision  frequency 


=  TTCTS  '-q 


and  f(v0)  =  Boltzmann  velocity  distribution 


=(*F )  v° 


exp 


-HVo' 


2kT 


with  a  =  hard  sphere  collision  diameter. 

Considering  the  case  of  small  energy  defects,  (14)  becomes 


(15) 


pr  n 
K*a  | 


S.E.D.L. 


16TT3n3LaU?BUa, 


(£-) 


naa  I(T) 


where  I (T)  is  the  following  integral 


®  /.  4LaEda  \ 

(16)  I  (T)  =  J  v*  y1 - ^“jexp 


uvo ' 


Ia(T)-Is  (T) 


2kT 


dv„ 


where 


00 

-^V0S 

J  Vo  eXp 

0 

2kT 

dvn 


l /-  5 


and 


I- (T)  =  4L2E 


aEf  J  v; 


o  exP 


uvr 


2kT 


dv. 


The  integrals  in  (16)  are  both  of  the  form 

1  • 3  *  5  *  *  *  ( 2n-l) 


(17) 


1 


x8  “  exp  l-ax* 


dx  = 


2n+1  a" 


TT  , 

a 


Therefore , 


(18) 


Ij  (T)  = 


m.  J 

2  11 


(*H 


and 


(19) 


Ia(T)(=  2/2  AaE2  (~“)3/4 


If  we  substitute  (18)  and  (19)  into  (15)  we  obtain  the  small 
energy  defect  rate  coefficient 


(20) 


R.r:l 

S.E.0.1. 


24/2  7 

=  -  "  "■  - 
h8 


nV2  L8  cr8(kT)%  U2n  U*B  |^1  - 


4uL2E? 
3kT 


This  rate  coefficient  can  be  simplified  into  convenient  units 
through  the  following  substitutions: 


(/-  6 


(  1.289-12) 


L3  (l”i6) 

of  (l-16) 

therefore. 


pr  ■ 

K.  i 


=  8.88"14 


'  S.E.O.L . 


of 


ip3 


U2 

r  n 


\  _  4hL3E2  -I 
3kT  I 


and 

=  5.703-  k*£ 

3kT  |_  T 


< 


where  Ea  is  in  cm”1 


The  small  energy  defect  rate  coefficient  then  becomes 


Rr  b 

«  a 


=  8.88”13  H, 


(22) 


S.E.D.L* 


Lf  of  T% 


«?.  Of.  - 


5.703”3 


LAEd 


The  large  energy  defect  rate  coefficient  can  be  similarly 
evaluated  as  follows.  If  (11)  and  (12)  is  inserted  into  (14)  we 
obtain 


(23) 


Kl 


IL.E.O.L. 


16ttVL3 

_ 


w 


(w) 


(T) 


where 


CD 

I0  (T)  =  J  v4  exp 


o 


-jiV3  _  4LEd 


2kT 


dv 


o  . 


(24) 


The  integral  in  (23)  cannot  be  evaluated  in  closed  form  but 
it  is  standard  for  this  type  of  integral  (A4)  to  expand  about 
the  maximum  of  the  exponent. 


The  maximum  of  the  exponent  is  obtained  as 


d_ 

dv. 


2kT 


4LE, 


=  0 


which  leads  to 


^4LEakT 


(25) 


Expanding  the  kernal  in  the  integral  l0 (T)  about  v„  we  obtain 

(v0-vH)8  aaKo 


Ko  (v0 )  =  K0  (vH  )  + 


+  H.O.T. 


The  second  derivations  in  (25)  can  be  evaluated  as 


(26) 


afi<o 

av8 


8LE, 
kT  Vg 


which  upon  substitution  of  vM  gives 
9aKo  .  -3 1 


av0a 


kT 


Substituting  (26)  and  (25)  into  (24)  yields 


(27) 


lo 


exp 


2kT 


4LEd 

vh 


exp 


(v0-v„)a 

~~2 


This  integral  is  now  of  a  standard  form  and  can  be  integrated  to 
give  (upon  extension  of  limits) 


(28) 


I0  (T)  =  v*  exp 


“liVj?  -  8kTLE4 
2kTvM 


4rTkT  \% 

—  J 
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The  exponential  function  in  (28)  can  be  simplified  via  vH  to 
give 


exp 

-^vj?  -  8kTLE4 

=  exp 

-6 

LEa 

2kTvM 

kT 

2 

Assembling  (29) , (28) ,  into  (23)  we  obtain 


which  upon  inserting  the  expression  for  v*  gives 


V 
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The  pre-exponential  factor  can  be  reduced  to  convenient  units 
via 


1 

6.02s3 


1.81-38 


(1“8) 


2. 14" 87 


a 


2 


) 


and 


-4/  4/ 

E4  /a  (2ttc)  /a 


1.091S 


and  where 


6.08-1* 


yielding 


=  1.508-1?T^  S  h/ 3  a*2  Ed  ^  USB 


And  the  exponent  can  be  reduced  similarly  to  yield 


V  -lO 


Thus,  the  exponent  becomes 


_  a/  \4 

-1.32  (LAEd)/a  j 

Inserting  the  pre-exponential  and  the  exponential  factor  in 
(30)  gives 


(31) 

7/  10/  _  4/ 

Kl  1  «  1.508-12  h/6l/3  a  a2  Ed  A 

3  Ur2n 

u.2. 

1  U  E*  D.l. 

—  ay 

3  /  M-H 

l  T 

y3 

•  exp  -1.32  (LAEd )  /: 

) 

Thus,  (22)  and  (31)  are  the  rate  coefficients  for  vibration- 
vibration  exchange  processes  in  the  small  energy  defect  and 
large  energy  defect  limits,  respectively.  The  units  for  these 
coefficients  can  be  listed  as  follows 


M-h 

La 

a  a 

T 

Ed 

and 


reduced  mass  in  molecular  units 

o 

exponential  parameter  in  A. 
hard  sphere  diameter  in 
temperature  in  °K 
energy  defect  in  cm"1 

matrix  element  for  r-«n  vibrational  transition. 
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FUNDAMENTAL  MOLECULAR  PROCESSES  IN  CO  LASER  SYSTEMS 


Introduction 

Despite  some  of  its  own  idiosyncrasies  (partial  inversion, 
multilevel  action)  the  CO  laser  has  appeared  to  be  a  tanta- 
lizingly  straight  forward  application  of  known  atomic  and  mo¬ 
lecular  processes  which  in  most  cases  have  been  theoretically 
and  experimentally  very  thoroughly  documented  in  recent  yea.rs . 
Indeed,  optimism  appears  to  be  reasonably  justified  in  view 
of  the  large  expertise  developed  for  analysis  of  more  complex 

reentry  and  upper  atmospheric  problems  which  involve  large 

1-3 

sets  of  relaxing  vibrational  modes  in  diatomics  .  .  Early 

4 

results  for  CO  and  CO/N2  systems  are  encouraging,  although 

experimental  results  involving  the  effects  of  minor  additives 

5-7 

such  as  Hg  addition  need  increased  study  for  inclusion  into 
models  for  molecular  laser  behavior. 

Fundamental  processes  worthy  of  consideration  are 'shown 
in  Figure  1.  Only  those  processes  which  create  the  laser  ac¬ 
tion  and  not  those  which  destroy  it,  are  included.  Other 
processes  which  must  be  considered  in  a  more  complete  treat¬ 
ment  of  this  system  would  include  ionization  of  CO  to  form 
C0+ ,  followed  by  dissociative  recombination: 

C0+  +  c  ->  C  +  0 

which  would  lead  tc>  carbon  deposits  and  presumably  eventual 
depletion  of  the  CO  over  long  term  operation.  Although  the 

\ 
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Figure  1 

FUNDAMENTAL  MOLECULAR  PROCESSES  IN  CO 
LASER  SYSTEMS 

1.  ELECTRON  IMPACT  EXCITATION 

e  +  CO  {CO“ }  •>  CO^  +  e 

* 

c  +  Hg  *►  Hg  +  e 


2.  ELECTRON  IMPACT  IONIZATION 
e  +  Hg  *►  Hg  +  2e 


3.  ELECTRONIC  -  VIBRATIONAL  ENERGY  TRANSFER 

*  4 

Hg  +  CO  -»  Hg  +  CO” 

4.  VIBRATION  VIBRATION  ENERGY  EXCHANGE 

N2(v)  +  CO(r)  -►  N2(v-1)  +  CO(r+l) 

5.  VI BRAT I ON -TRANSLATION  ENERGY  TRANSFER 

CO(r)  +  He  -►  CO(r-l)  +  He 

6.  RADIATIVE  PROCESSES 

Spontaneous 

Induced 
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V 


rate  constants  for  this  particular  reaction  are  not  at  present 

+  +  4. 

known,  analogous  reactions  for  N'O  »  ^  ,  etc.  are  well 

g 

established  .  The  cross  section  data  for  electron  impact  ex- 

9  10 

citation  of  both  diatoinics  and  atoms  has  been  available  for 
some  time  and  recently,  the  rate  constants  for  such  processes 
in  N2  and  CO  have  been  documented^. 

In  the  case  of  Hg,  or  any  low  ionizaiton  potential  addi¬ 
tive,  an  important  result  of  the  ionization  step  could  be  addi¬ 
tional  electrons  with  a  resulting  electron  temperature  decrease 
due  to  the  thennalization  of  the  hot  electrons  by  the  near  ther¬ 
mal  electrons  produced  in  the  Hg  ionization.  If  one  examines  the 
limiting  case  of  Iig+  equilibrium  with  the  electron  temperature 
through  the  Saha  equation,  it  is  highly  probable  at  T  =  20,000°K 
that  most  of  the  Hg  will  be  ionized.  At  a  vapor  pressure  of  Hg, 

1  TJ 

of  lp,  this  could  represent  10  "  elcctrons/cc  as  an  upper  limit. 
The  remaining  Hg  (or  Xe)  which  is  not  ionized  will  be  excited 
electronically,  and  despite  the  fact  that  the  radiative  lifetime 
of  Hg  (  P)  is  10  sec,  radiation  trapping  will  cause  substantial 
quenching  by  CO  or  N2  with  resultant  conversion  of  electronic 
energy  to  vibrational  energy. 


Electron  Impact  Processes 

Consider  first  the  excitation  of  by  electron  impact. 
Because  this  process  proceeds  through  an  unstable  negative  ion 


intermediate  state,  then 


activation  energy  for  the  pro¬ 


cess  which  leads  to  a  threshold  of  about  1.5  ev  rather  than  the 


0.29  ev  vibrational  spacing.  Cross  sections  are  large  for  col¬ 
lision  energies  between  2  and  4  ev,  approaching  gas  kinetic  at 
.the  maximum,  for  excitation  from  the  ground  vibrational  level 
to  levels  from  1  through  8  as  shown  in  Figure  2.  Beyond  6  or  7, 
the  cross  sections  diminish  rapidly  and  are  therefore  of  little 
importance  compared  to  excitation  from  excited  vibrational  levels 
to  still  higher  levels  as  will  be  seen. 

Similar  data  exist  for  CO  as  shown  in  Figure  3.  It  is  ' 
significant  that  the  threshold  for  CO  excitation  occurs  below 
1  ev  and  that  the  cross  sections  are  four  times  as  large  for 
CO  as  for  N2  excitation. 

Aside  from  putting  energy  into  the  vibrational  mode,  the 
implications  of  these  electron  impact  processes  are  illustrated 
in  Figure  4,  from  Nighan's  work  .  The  electron  energy  (u)  dis¬ 
tribution  f  (u)  =  2  J=-  Te  e  U^Te  would  be  a  straight  line  if 

the  distribution  were  Maxwellian.  The  large  deviations  from  a 

% 

straight  line,  characterized  by  the  absence  of  electrons  with 
energy  greater  than  about  2  ev  for  N2  is  the  result  of  the  large 
vibrational  excitation  cross  sections  which  provides  an  energy 
loss  mechanism  nearly  preventing  electrons  from  reaching  energies 
above  2  ev. 

However,  as  Ty  increases  the  reverse  mechanism  rebuilds 
part  of  the  distribution,  and  this  tail  may  in  fact  be  important 
for  electronic  and  ionization  processes. 

However,  because  we  are  interested  in  rate  constants  or 
integrated  cross  sections,  we  need  not  worry  about  deviations 
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F  i  g  .  2  Energy  de¬ 
pendence  of  the  Vibra¬ 
tional  cross  section,  of 
nitrogen  by  electron  im¬ 
pact.  The  curves  are 
obtained  from  sets  of 
curves  similar  to  Fig.  3 
Utica  at  different  ener¬ 
gies  of  the  incident  elec¬ 
trons.  When  the  ordi¬ 
nate  numbers  are  mul¬ 
tiplied  by  10~S*  a  cross- 
section  scale  (in  cm1)  is 
obtained  such  as  to  give 
a  total  vibrational  cross 
section  in  agreement 
rith  Haas.  See  text  for 
a  discussion  of  errors  in 
aoss-section  scale. 


calculated  for  various  vibrational  temperatures 


e-  6 


from  a  Boltzmann  distribution  for  vibrational  excitation  as 
illustrated  in  Figure  5.  Rate  constants  for  ^  excitation 
are  insensitive  to  the  electron  kinetic  temperature  above  1 
ev  or  11,000°K.  At  this  point,  the  rate  constants  are  very 
large.  It  is  especially  interesting  to  note  the  comparison 
between  kQ  4  and  k^  g  which  at  high  Tg  are  very  n’early  equal. 
Thus,  at  high  vibrational  temperatures  a  two  step  process  may 
be  the  preferred  electron  excitation  route.  Again,  the  tail 
of  the  electron  distribution  is  much  more  important  in  exciting 
electronic  states  of  additives  than  for  increasing  the  vibra¬ 
tional  excitation  rates. 

Quenching  of  Electronic  States  by  Molecules 

Once  an  additive  atomic  species  has  been  electronically 
excited,  there  is  the  possibilities  of  energy  transfer  to  the 
vibrational  mode  depending  upon  the  potential  energy  curves  of. 
the  system  quenching  system.  Two  types  of  quenching  mechanisms 
are  shown  in  Figure  6.  First  is  the  ionic  intermediate  mecha¬ 
nism,  which  is  the  preferred  route  for  the  quenching  of  excited 
electronic  states  of  low  (<  4  ev)  ionization  potential  atomic 
species  such  as  the  alkali  atoms,  by  diatomic  molecules,  two 
important  examples  of  which  are  ^  and  CO.  The  second  is  a  co¬ 
valent  mechanism  which  proceeds  directly  through  a  curve  crossing 
mechanism  involving  the  direct  crossing  of  reactant  excited  state 
of  a  stable  or  mctastable  intermediate  and  product  potential 
curves.  An  example  of  this  type  is  the  quenching  of  0(*D)  by  ^ 
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Figure  6 


QUENCHING  OF  ELECTRONIC  STATES  BY  MOLECULES 

* 

A  +  M2  A  +  M2 

Ionic  Curve  Crossings 

Na(32P)  +  N2  -  {Na+  -  N2‘}  ->  Na(32S)  +  N 

Covalent  Curve  Crossing 

0(1D)  +  N2  -»•  {N20}*  ->  0(3P)  +  N2 

Hg(63Pij0)  +  CO  -  {HgCO}*  -  HgC61S0)  +  C 


The  principle  involved  in  the  ionic  mechanism  is  illus¬ 
trated  for  atom-atom  quenching  in  Figure  7.  At  large  interac¬ 
tion  distances,  and  hence  with  potentially  large  cross  sections 
the  potential  curves  corresponding  to  the  ionic  intermediate 
complex  become  energetically  favorable  compared  to  the  excited 
curve  providing  a  curve  crossing  mechanism  coupling  initial 
(reactant)  curves  with  final  (product)  curves.  The  curve  cross¬ 
ing  at  each  point  of  intersection  can  be  independently  evaluated 
and  is  assumed  to  be  given  by  the  Landau  Zener  form  .  The  tran¬ 
sition  matrix  element  is  obtained  from  the  electronic  component 
(known  from  the  Hasted  Chong  correlation)*^.  The  vibrational 
overlap  integral  at  the  crossing  point,  which  is  obtained  from 
the  appropriate  Franck  Condon  factors,  is  necessary  to  describe 
excited  atom  quenching  by  molecules.  The  potential  energy  curves 
for  electronic  states  have  superimposed  on  them  the  vibrational 
levels  of  both  the  diatomic  and  the  negative  ion,  so  that  the 
overall  partial  cross  sections  are  computed  via  a  classical  dif¬ 
fusion  through  the  illustrated  potential  grid  (Figure  8). 

2 

Sample  results  of  Na  (5  P)  quenching  by  are  shown  in 
Figure  9.  The  significant  feature  of  the  result  is  the  non-reso- 
nant  behavior  and  the  large  fractional  energy  which  goes  into 
vibrational  (over  SO') . 

The  potential  energy  curves  for  a  covalent  quenching 
reaction  such  as  0(^D)  +  ^  are  sketched  in  Figure  10*^.  For 
this  case,  the  interaction  matrix  element  at  each  crossing  point 
is  given  by  an  electronic  contribution  due  to  spin-orbit  coupling 
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FIGURE  *1.  Potential  Energy  Curves  for  an  Atom-Atom  System 
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Figure  8:  Simplified  potential  curves  for  an  excited  atom 
(  molecule  quenching  reaction.  i 
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FINAL  SURFACES 


while  the  vibrational  overlap  is  to  compute,  to  lowest  order, 
between  vibrational  wave  functions  of  the  same  molecule,  i.e. 

■  N2  in  this  case.  Thus,  the  only  quenching  channels  which  have 
significant  interaction  elements  (o.r  cross  sections)  are  those 
in  which  the  vibrational  state  of  the  reactant  N2  is  nearly  un¬ 
changed.  Thus,  the  expected  vibrational  excita-tion  produced  in 
covalent  quenching  reactions  is  small*^.  This  result  has  been 

indirectly  verified  by  recent  atmospheric  vibrational  temperature 

.  .  16 
analysis 

In  the  case  of  Hg  quenching  by  CO,  the  potential  energy 
curves  arc  not  available  for  theoretical  analysis,  but  experi¬ 
mental  data  is  available.  The  system  is  shown  in  Figure  11,  and 

1 7 

the  experimental  results  of  Polanyi  arc  shown  in  Figure  12. 

*  -  7 

Even  though  Hg  lifetimes  are  10  sec,  as  shoitfn  in  Figure  13,  it 
is  well  known  from  very  early  studies  on  Hg  areas  that  significant 
radiation  trapping  occurs  under  plasma  conditions  so  that  an  im¬ 
portant  fraction  of  the  electronic  energy  eventually  channels 
into  the  vibrational  mode. 

Vibrational  Exchange 

The  vibrational  Master  equation*  is  illustrated  in  Figure 

14.  A  slowly  time  varying  quasistcady  state  solution  to  the 

1 8 

Master  has  been  found  by  Trcanor  to  be 

nr  =  no  exp  ^"rY^  exp  C'Er/kT) 

=  nQ  exp  (-Er/k°r*), 
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Partial  Quenching  Cross  Sections  in  Ang. 
for  0(*D)  +  N2  -  co-linear  case 


triplet  vibr. 
level 

Maximum  Spin -Orb it 

300  cm'1 

Interact 

500  cm'1 

6 

0.0 

‘0.0 

5 

0.0 

0.0 

4 

0.0 

0.0 

3 

0.0 

0.0 

2 

0.0 

0.002 

1 

0.0 

0.33 

0 

2.8 

5.7 

total  quenching  cross 

section 

2.8 

6.03 

B~  16 


LOg  \'y  ( ARBITRARY  VNTTS5 


•  <000  eooo  12030  i:033  JOOOO 


Ev  (cm*') 


Viliration.il 
distribution  of  CO, 
showing  the  effect 
of  impurity  in  CO. 
(a)  Run  (1) : impure 
CO  without  diluent 
+,  (1>)  1-tun  (0) : 
impure  CO  with  Ar 
O,  (c)  I  tun  (2): 
purif.ed  CO  withuut 
diluent  O. 


Fig.  12a 


Vibrational 
distribution  of  (pur¬ 
ified)  CO,  showing 
the  effect  of  inert 
diluent  with  li i”li  and 
with  normal  mer¬ 
cury  pressure,  and 
the  effect  of  ‘vibra¬ 


tional  scavenger’,  (a) 
Kun  (>!):  CO  in  Ar, 
high  Jig  pressure, 
□.  (b)  Run  (5): 
CO  in  Ar,  normal 
llg  pressure,  O,  (c) 
Run  (3) :  CO  in  Nj, 


normal  Hg  pressure, 
A,  (d)  Run  (2):pu- 
liftcd  CO  without  dil¬ 


uent,  dashed  cirelc. 
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Fig.  12b 
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Figure  14 

RELAXATION  BY  VIBRATION  EXCHANGE. 


master  equation 

"Ati,  [»r+i»-  exp  [»,»-  exp  «,-ikJ  +  X)  Pf+iA1 

(  F,»+E,-1-Er-E,\  1  V' n  /  E,-\E.-Enl~Er^ 

X  J^K, +!«,-!-  exp  ^ - — - J  »»,«,!  -  X)  ■Pf.r-l*’*1  I  exp  y - — - 

+0f+l.f  [^Mr+lA7-  exp  «rtfj  [»,#-  exp  — ■ j~- j 


+  £<W-“  «<p(-  ^±^P^)m'»] 

-  E  &/-/■*”  £*V>’*-  e«P  (-  ^+F/i  ty~‘  j  • 
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)!l+xn,-.i 


Nr  *=  Nq  exp  (-Ry)  exp  (-FR/kT) 

=  Nq  exp  C-rR/k4>R*)  , 

where  it  may  be  shown  by  substitution  that  y  is  the  same  for 

*  -  * 

both  molecules,  and  where  0r  and  <J>R  are  population  factors 
for  the  rth  and  Rth  vibrational  levels  of  species  A  and  B, 
respectively,  y  may  be  evaluated  in  terms  of  level  1  of  either 
molecule : 

y  =  (Kj/k)  Cl/e^-T'1)  =  (Fj/k)  (l/^-T'1), 


so  that 


-  E  rE  0,  rF. 

n  =  n  exp  { - -J  [-g-  -  -y-  (nr"  '  !)]> 

k0j  r  r 


and 


Nr  =  No  exp  {-  [ 


RF. 


V- 


k<}>. 


T 


R 


RF 

C-~  -  1)]} 

FR 


For  anharmonic  oscillators,  these  distributions  are  clearly 
non-Boltzmann  since  the  terms  in  brackets  are  dependent  on  the 
vibrational  level.  In  the  harmonic  model,  Ef  =  rEg  and  FR  = 
RFj ,  therefore 


nr  =  no  exp  (-rEg/kOjj)  , 


and 


Nr  =  No  exp  (-RFj/k^j)  , 
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where  0^  =  0^  and  =  4*2  »  are  the  vibrational  temperatures 
for  species  A  and  B,  respectively.  Here,  we  find  Boltzmann 
distributions  with  the  relationship  between  0^  and  4^  given  by 

(E1/eH)  -  (F^n)  -  CEX 

Further,  the  relationship  between  populations  in  the  same  level 
R  is 

nR  n  Z  R(E,  -  F,) 

NJ  *  2  exp  [ - YT~ - 1  » 

where  n  and  N  are  the  total  concentrations  and  z  and  Z  are  the 
vibrational  partition  functions  of  species  A  and  B,  respectively, 
defined  by 

z  ~  E  exp  (-  REj/kOjP 
R 

and 

Z  *  E  exp  C-RF1/k<f.H) . 

R 

It  is  important  to  note  that  the  VV  steady-state  forces  the 
relative  distributions  of  the  two  species  to  be  nearly  charac¬ 
terized  by  the  kinetic  temperature  even  though  the  vibrational  ■ 
temperature  of  each  species  may  be  considerably  higher  than  given 
by  equilibrium.  The  ratio  of  partition  function  for  species  A 
and  B  is  weakly  dependent  on  the  vibrational  temperature,  and 
thus,  deviates  slightly  from  equilibrium.  This  characteristic 
of  steady-state  VV  exchange  processes  in  gas  mixtures  in  the  har¬ 
monic  model,  is  also  found  in  the  anharmonic  model.  This  is  the 


important  pumping  mechanism  which  creates  the  partial  inversion 
in  the  CO  laser. 
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In  Figure  15a  we  show  the  LeGay  demonstration  of  the 
validity  of  the  VV  model  for  N,  -  CO  systems.  In  Figure  15b, 
the  effect  of  rotational  temperature  which  will  make  the  agree¬ 
ment  of  Figure  15a  better  is  shown.  Other  verifications  have 

20  21  22 
been  done  by  Cornell  ,  Aerospace  ,  (Figure  16),  AVCO  ,  etc. 

A  time  dependent  model  for  a  pulsed  CO  system  including 

all  these  energy  transfer  processes  for  a  pure  CO  system  of  20 

levels  is  shown  in  Figure  17.  The  vibrational  population  factor 

or  temperature  for  each  vibrational  level  in  a  20  level  program 

in  CO  and  He  is  shown  as  a  function  of  time,  as  computed  through 

1  3 

the  Runge  Kutt.a  Merson  *  integration  technique  including  elec¬ 
tron  impact  excitation  for  all  levels,  VV,  VT  processes  and  spon¬ 
taneous  emission.  A  comparison  with  the  Rich^  steady  state 
calculation  is  shown. 

« 

A  similar  model  describing  the  pulsed  behavior  of  ^  -  CO  - 
He  in  the  presence  of  hot  electrons  has  also  been  developed. 

The  presence  of  ^  plays  only  a  small  role  in  enhancing  the  CO 
vibrational  populations  unless  the  ^  concentration  is  greater 
than  the  CO. 
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F  i  p  .  15  a.  Rearranged  Boltzmann  te'mpcrature  vs.  CO  relative  concentration :  curves  given  by  Fisher  and  Kummfer 

(1968b)  for  r,°(Nj)  =  2000  °K.  kinetic  temperature: - ,  T  =  300  “lC; - ,T  -  500  °K. - -.curve  given  by 

Fisher  (private  communication)  for  TP°( Nj)  *=  2200  °K,  T,c(CO)  ~  300  °K,  T  —  400  °K.  Oi  present  experimental 
points  (r,°(Na)  ~  2200  °K)  T  =  400  °K. 


Fig.  1  5  b  .  Vibrational  energy  Inst  inio  rotational  and  kincti 
energies  from  initial  vibraiional  distribution  to  steady’  state  du 
to  ronresonant  vibrational  exchange  processes  as  a  function  c 
composition  for  X,  and  CO  mixtures  at  different  initial  X,  vibra 
uonal  temperatures  and  a  kinetic  temperature  of  300’K  Th 
points  shown  are  calculated. 
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Appendix  III 


THE  VIBRATIONAL  EXCITATION  OF  CO  BY 
COLLISION  WITH  HELIUM  * 

R,  MARRIOTT 

Research  Institute  for  Engineering  Sciences 
Wayne  State  University,  Detroit,  Michigan,  U.S.A. 

Inelastic  processes  involving  the  vibrational  excitation 
of  molecules  by  collision  with  inert  gases  are  of  particular 
relevance  to  the  study  of  gas  laser  phenomena.  In  this  paper 
cross  sections  for  the  vibrational  excitation  of  CO  by  Ke 
impact  are  calculated,  using  a  computer  code  previously 
applied  to  the  case  of  CO-CO  collisions (1) . 

The  code  is  based  upon  a  close  coupling  analysis  of  the 
collision  system.  Rotational  states  of  the  target  and 
internal  structure  of  the  impacting  particle  are  neglected. 
The  molecular  wave  function  is  assumed  to  be  separable  into 
electronic,  rotational  and  vibrational  components.  The 
scattering  potential  is  similarly  assumed  to  be?  separable 
in  terms  of  the  vibrational  and  intermolccul ar  coordinates 
and  to  be  spherically  symmetric.  An  empirical  potential 
function  is  subsequently  substituted  for  the  intcrmolecular 
component . 

To  this  approximation  the  radial  partial  vrrvc  equations 
are  solved  numerically  with  coupling  rat-'  a.  ’  between  any 
specified  number  of  molecular  states.  Total  cross  sections 
are  obtained  by  graphical  integration  over  all  significant 
partial  waves. 

Previous  calculations  have  found  the  Lennard- Jones 
function  4e  ( (cr/r)  2-  (o/r)  to  give  a  sal:  f  letory  repres¬ 
entation  of  the  scattering  potential.  T!v .•  n„:  ■ '.metric  values 
required  for  the  CO-IIe  system  have  b^c.n  tabulated  by 
Hirscbfelder  ct  al.(2):  e  =  9 . 5xl0~ J  (a . u . )  ;  0  =  5.93(a^). 

The  remaining  parameters  required  as  data  for  the  code 
are  the  oscillator  reduced  mass  (1.2Cxl04n.u.) ,  the  level 
separation  for  CO  vibrational  states  (0.2GG  eV)  and  the 
coll  ision  reduced  mass  (G  .  4  3x10 -’a  .  u . )  . 


Total  cross  sections  for  the  excitation  o!  CO  from  t he 
ground  to  th.  first  vibrational  state  resulting  from  fchs 
use  of  the  I.cnnard-Joncr  potential  (LJ)  are  shewn  in 
tabic  1.  Coupling  between. up  to  four  vibrational  states 
was  retained.  An  indication  of  the  r.vno.-j  o*J  acoi  racy  is 
given  by  t  he  rocy 

’  1  l  t’lTe  n  i  gtii  ficant.  fj cures. 
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.  While  there  are  no  direct  measurements  of  these  cross 
sections  for  comparison  with  the  calculated  values,  they 
can  be  used(l)  to  determine  the  vibrational  relaxation  time 
of  CO  in  He,  for  which  experimental  data  is  available (3) . 

Assuming  only  nearest  neighbor  transitions  to  be  import¬ 
ant  the  relaxation  time  can  be  written:  t  =  1/(y-j  sec- 

where  the  required  rate  coefficient  is  given  by 

Y01’=  4(10)11{MT)"J5|q01(E)  (E/kT)e“E/kTd(E/kT) 

M.  a.u.  being  the  collision  reduced  mass,  T  °K  the  gas  temp¬ 
erature  and  Qqi (E)  is  the  excitation  cross  section  for  the 
0-*l  transition  at  a  collision  energy  of  E  in  units  of  (Tra^)  . 

0 

The  theoretical  relaxation  times  obtained  in  this  way 
are  shown  in  figure  1  in  comparison  with  the  experimental 
values  of  Hooker  and  Millikan(3).  It  is  immediately  apparent 
that  theory  is  consistently  less  than  experiment  by  a  factor 
of  fifty,  indicating  that  the  cross  sections  are  too  large 
by  the  same  amount.  This  in  turn  implies  that  the  Lennard- 
Jones  repulsive  core,  which  dominates  the  excitation  process, 
is  too  strong. 

The  scattering  studies  by  Amdur  et  al.(4,5)  of  a  number 
of  gas  molecule  potentials  indicate  that  the  repulsive  part 
of  the  interaction  between  CO  and  the  inert  gases  behaves 
as  r~7  rather  than  r”1  . 

Although  the  CO-He  interaction  has  not  been  measured 
directly,  data  is  given  on  the  ^-He  system  and  it  is  shown 
that  the  CO  and  interactions  are  very  similar. 

On  this  basis  we  estimate  the  repulsive  interaction  for 
CO-He  to  be  232/r'(a.u.)  for  separations  3  <  r(ap)  <  4.3. 

For  the  long  range  part  the  Lennard-Jones  potential  is 
assumed  still  valid  and  the  two  functions  are  smoothly 
joined  over  a  small  region  by  the  common  tangent: 

V (r)  =  .04739-. 009158r  for  4.53  <  r(aQ)  <  4.84. 

Use  of  this  soft  core  (SC)  potential  reduces  the  strength 
of  the  coupling  between  vibrational  states  by  about  a  factor 
of  two.  Representative  values  of  the  SC  cross  sections  are 
shown  in  table  1.  Coupling  has  been  retained  between  up  to 
six  vibrational  states. 

Recalculation  of  the  vibrational  relaxation  time  on  the 
basis  of  these  crocs  sections  leads  to  the  third  curve  shown 
in  figure  1.  The  theoretical  results,  are  now  within  a  factor 
of  three  of  experiment,  except  at  temperatures  be  lev.*  300  UK 
where  there  is  a  sharp  decrease  in  the  temperature  dependence 
of  the  calculated  relaxation  time. 


Table  1 

Cross  sections  for  the  vibrational  excitation  of  CO  by  He. 


Potential 

LJ 

SC 

SC 

SC 

Collision  * 
energy  (cV) 

Excitation  cross 

Q01  Q01 

sections 

°02 

g12 

0.266 

0.600 

1.000 

1.500 
2.000 

2.500  . 

0.0** 

i.6"f 

2.4": 

9.7"1 

1.9 

2.6 

0.00  . 
1.00, 
1.28": 
1.13": 

4.06", 

8.72"-1 

_  o 

0.77  l 

2.60"J 

1.66, 

1.89", 

8.58"J 

6.72“® 

4.40", 

2.82", 

8.68"f 

1.64"-1 

*  Collision  energy  relative  to  the  ground  state. 
**  Superscript  indicates  power  of  10  by  which 
number  must  be  multiplied. 


Fig.  1  Vibrational  relaxation  of  CO  in  He  at  one  atmosphere. 
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